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13  ABSTRACT 


The  recently  observed  exponential  frequency  dependence  of  the  infrared  absorption  coeffi¬ 
cient  P  in  NaCl-  and  CaF^- structure  crystals  is  explained  in  terms  of  multiphonon  absorption 
processes.  A  calculation  gives  the  correct  frequency  dependence  and  approximate  magnitude 
of  P  and  predicts  the  temperature  dependence,  which  should  be  observable.  Use  of  the  central- 
limit  theorem  reduces  the  mathematical  complexity  of  the  analysis  sufficiently  to  yield  a  simpl^ 
closed-form  approximation  to  p.  In  the  exponential  region,  the  difference  processes,  in  which 
some  thermally  excited  phonons  are  annihilated,  are  negligible  with  respect  to  the  summation 
processes.  Terms  higher  than  second  order  in  the  perturbation  expansion  are  found  to  be  neg¬ 
ligible  for  small  n  (small  number  of  phonons  created).  An  explanation  involving  finite  phonon 
lifetimes  is  proposed  to  explain  the  fact  that  the  alkali  halides  show  less  structure  in  the  p-oc 
curves  than  do  the  semiconductor  crystals.  Preliminary  results  on  a  new  approach,  using  an 
anharmonic  Einstein  model,  also  exhibit  the  exponential  frequency  behavior..  In  addition, 
p~  cx)'2  is  predicted  for  "square-well"  type  impurities. 

Small  amounts  of  macroscopic  inclusions  in  or  on  the  surface  of  a  nonabsorbing  host  crys¬ 
tal,  in  some  cases  a  fractional  volume  of  ~10"7  to  10"8,  can  give  rise  to  a  value  of  the  optical 
absorption  coefficient  p  of  lO'^cm'l  (a  typical  value  of  current  interest).  For  various  types 
of  inclusions,  the  frequency  dependence  ranges  from  increasing  as  up,  to  independent  of  to,  to 
exponentially  decreasing  with  U).  The  temperature  dependence  ranges  from  independent  of  T 
to  increasing  as  TP  in  the  high-temperature  limit,  where  p  -  2  -  4  typically.  Simple  expres¬ 
sions  for  the  absorption  cross  section  are  derived  for  various  cases  of  practical  interest.  The 
cross  sections  are  used  to  derive  expressions  for  P  for  the  four  cases  of  large  inclusions  of 
strong  and  weak  absorbers  and  of  small  inclusions  of  dielectric  and  metallic  particles. 
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5  A  result  of  general  interest  in  the  theory  of  scatter¬ 
ing  is  that  the  absorption  cross  section  of  large  spheres 
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energy  densities  as  low  as  ~5  J/cm^  (excluding  Bloem- 
bergen  field  enhancement)  can  cause  material  failure. 

Experimental  curves  of  ^  vs  uj,  collected  from  the 
literature,  are  presented  for  a  number  of  materials. 
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sivity  measurements,  effects  of  pressure  on  the  opera¬ 
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sion  of  the  previously  observed  T^  linewidth  dependence 
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technical  progress  reports,  are  summarized  briefly  in 
the  Introduction. 
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PREFACE 


This  Final  Technical  Report  describes  the  work  performed  on  Contract 
DAHC15-72-C-0129  on  Theoretical  Studies  of  High-Power  Infrared  Window 
Materials  during  the  period  from  July  1,  1972  to  November  30,  1972.  Two 
previous  technical  progress  reports  were  prepared  under  this  contract. 

Other  publications  are  listed  in  the  Introduction  (Sec.  A). 

The  work  on  Theoretical  Studies  of  High-Power  Infrared  Window  Materials 
will  continue  under  Contract  DAHC15-73-C-0127. 

Computer  analyses,  which  have  been  delayed  by  the  change  from  6600  to 
7600  computers  at  the  laboratory  of  Dr.  Arnold  Karo  (consultant  from  Liver¬ 
more)  are  now  under  way.  The  results  will  be  presented  in  a  future  report. 

The  following  consultants  are  participating  in  the  program  at  present: 

Dr.  Arnold  M.  Karo,  Lawrence  Livermore  Laboratory,  Dr.  Alexei  A. 
Maradudin,  University  of  California,  Irvine,  Dr.  Douglas  L.  Mills,  University 
of  California,  Irvine,  and  Dr.  Lu  J.  Sham,  University  of  California,  San  Diego. 


SUMMARY 


The  recently  observed  exponential  frequency  dependence  of  the  imrared 
absorption  coefficient  j3  in  NaCl-  and  CaF2_ structure  crystals  is  explained  in 
terms  of  multiphonon  absorption  processes.  A  calculation  gives  the  correct 
frequency  dependence  and  approximate  magnitude  of  0  and  predicts  the  tem¬ 
perature  dependence,  which  should  be  observable.  Use  of  the  central-limit 
theorem  reduces  the  mathematical  complexity  of  the  analysis  sufficiently  to 
yield  a  simple  closed- form  approximation  to  0.  In  the  exponential  region, 
the  difference  processes,  in  which  some  thermally  excited  phonons  are  anni¬ 
hilated,  are  negligible  with  respect  to  the  summation  processes.  Terms 
higher  than  second  order  in  the  perturbation  expansion  are  found  to  be  neg¬ 
ligible  for  small  n  (small  number  of  phonons  created).  Ar.  explanation  in¬ 
volving  finite  phonon  lifetimes  is  proposed  to  explain  the  fact  that  the  alkali 
halides  show  less  structure  in  the  (3 -to  curves  than  do  the  semiconductor 
crystals.  Preliminary  results  on  a  new  approach,  using  an  anharmonic 

Einstein  model,  also  exhibit  the  exponential  frequency  behavior.  In  addition, 
-2 

0  ~  w  is  predicted  for  "square-well"  type  impurities. 

Small  amounts  of  macroscopic  inclusions  in  or  on  the  surface  of  a  non- 

-7  "8 

absorbing  host  crystal,  in  some  cases  a  fractional  volume  of  ~10  to  10  , 

-4  -1 

can  give  rise  to  a  value  of  the  optical  absorption  coefficient  j3  of  10  cm 

(a  typical  value  of  current  interest).  For  various  types  of  inclusions,  the  fre- 

2 

quency  dependence  ranges  from  increasing  as  to  ,  to  independent  of  to,  to 
exponentially  decreasing  with  co.  The  temperature  dependence  ranges  from 


1 


independent  of  T  to  increasing  as  T*3  in  the  high-temperature  limit,  where 
p  s  2  -  4  typically.  Simple  expressions  for  the  absorption  cross  section  are 
derived  for  various  cases  of  practical  interest.  The  cross  sections  are  used 
to  derive  expressions  for  j3  for  the  four  cases  of  large  inclusions  of  strong 
and  weak  absorbers  and  of  small  inclusions  of  dielectric  and  metallic  particles. 

A  result  of  general  interest  in  the  theory  of  scattering  is  that  the  absorp- 

2 

tion  cross  section  of  large  spheres  is  equal  to  ( 1  -  R  •  )tra  ,  where  Rj^  is  the 

2 

reflection  coefficient  for  normal  incidence,  rather  than  (1  -  (R))7Ta  ,  where 

(R)  is  the  geometric -optics  average  reflection  over  the  illuminated  hemisphere. 

Material  failure  resulting  from  local  heating  of  macroscopic  inclusions  is  a  far 

greater  problem  in  high -intensity,  short-pulse  systems  than  in  low-intensity 

long-pulse  or  cw  systems  having  equal  average  intensity.  Microsecond  pulses 

2 

with  energy  densities  as  low  as  ~5J  /cm  (excluding  Bloembergen  field  enhance¬ 
ment)  can  cause  material  failure. 

Experimental  curves  of  0  vs  u>,  collected  from  the  literature,  are  pre¬ 
sented  for  a  number  of  materials.  Results  on  multiphonon  calculations,  the 
need  for  emissivity  measurements,  effects  of  pressure  on  the  operation  of 

infrared  windows,  effects  of  interband  and  deep-level  transitions  on  infrared 

2 

absorption,  and  a  discussion  of  the  previously  observed  T  linewidth  dependence 
for  NaCl,  which  were  presented  in  the  two  previous  technical  progress  reports, 
are  summarized  briefly  in  the  Introduction. 


2 


Sec.  A 


A.  INTRODUCTION 


In  the  first  technical  progress  report  the  program  was  outlined,  the 
problems  and  objectives  were  defined,  and  the  relevance  to  other  Department 
of  Defense  programs  was  discussed.  A  highly  simplified  discussion  of  the 
theory  of  infrared  absorption  was  given  as  an  introduction  for  investigators  in 
1  elated  fields.  Preliminary  calculations  of  the  values  of  the  absorption  coef¬ 
ficient  /3  for  the  two-phonon  and  the  general  n-phonon  summation  intrinsic 
processes  provided  an  overview  of  the  theoretical  problem  and  established  the 
method  of  calculation  and  the  types  of  approximations  necessary  to  render  the 
mathematics  tractable.  Methods  of  estimating  the  values  of  the  anharmonic 
coefficients,  the  higher-order  terms  in  the  dipole-moment  expansion,  and 
the  joint  densities  of  states  were  established.  The  exponential  decay  of  /J  as 
a  function  of  frequency  was  explained  in  terms  of  multiphonon  absorption  pro¬ 
cesses.  It  was  suggested  that  two -phonon  summation  processes  in  materials 
having  the  NaCl  structure  are  small  unless  one  phonon  is  optical  and  the  other 

is  acoustical.  The  latter  result  afforded  a  tentative  explanation  of  the  previously 
2 

observed  T  dependence  of  the  resonance  linewidth  of  NaCl  and  the  small 
number  of  two -phonon  summation  peaks  in  various  alkali  halides. 


It  was  pointed  out  that  currently  known  transmission  and  calorimetric 
measurements  cannot  be  used  to  obtain  values  of  jS  in  a  crucial  range  of  fre¬ 
quencies,  and  that  emissivity  measurements  should  be  capable  of  providing  the 
needed  values  of  p  in  the  full  range  of  interest.  Pursuit  of  such  experimental 
programs  would  contribute  substantially  *o  the  information  presently  needed  for 
Department  of  Defense  laser  programs.  Two  effects  of  pressure  on  the  operation 


3 


Sec.  A 


of  windows  were  considered.  First,  the  deformation  resulting  from  a  static 
pressure  differential  between  the  two  faces  of  the  window  can  distort  the  laser 
beam  or  fracture  the  window.  Second,  the  window-material  figures  of  merit 
for  thermally  induced  optical  distortion  that  were  previously  developed  are 
changed  when  the  thickness  required  to  withstand  the  pressure  is  taken  into 
consideration. 

It  was  pointed  out  that  several  references  on  inter-valence-band  transi¬ 
tions  and  transitions  involving  deep  levels  in  the  band  gap  for  semiconductor 
materials  suggested  that  at  impurity  concentrations  so  low  that  the  usual  free- 
carrier  absorption  is  not  detectable,  transitions  of  holes  between  two  different 
valence  bands  apparently  can  give  rise  to  broad  absorption  bands  that  can  be 
strong  (/3  =  0.3  cm  in  one  case,  for  example)  in  the  2-6fim  region.  Transi¬ 
tions  of  electrons  between  different  regions  of  the  conduction  band  where  the 
densities  of  states  are  high  also  should  give  rise  to  observable  values  of  /3. 

3 

In  the  second  technical  progress  report,  the  results  of  a  Green's  function 
analysis  and  additional  ordinary  perturbation -theory  analyses  of  the  n -phonon 
summation  process  were  presented.  A  detailed  analysis  of  the  multiphonon 
absorption  was  presented  and  applied  to  several  alkali  halides.  It  was  sug¬ 
gested  that  a  possible  explanation  of  the  fact  that  the  j3  -  to  curves  of  the  alkali 
halides  and  alkaline -earth  fluorides  show  less  structure  than  those  of  the  semi¬ 
conductor  materials  is  that  the  greater  anharmonicity  of  the  NaCl- structure 
crystals  could  give  rise  to  such  short  lifetimes  of  the  zone-boundary  phonons 
that  the  peaks  are  broadened  to  the  point  of  essentially  complete  overlap. 
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The  contributions  to  the  n -phonon  absorption  coefficient  ^  from  higher- 
order  processes  with  n  phonons  in  the  final  state  were  calculated  for  a  number 
of  processes  and  classes  of  processes.  Terms  containing  self-energy  diagrams 
within  the  final-state-  or  fundamental-phonon  lines  simply  renormalize  the 
phonons.  Diagrams  having  one  intermediate  phonon  and  one  3-phonon  vertex 
give  the  leading  contributions  to  of  the  higher-order  processes  considered. 

In  NaCl  at  300  K,  the  contributions  to  <3^  from  the  higher-order  process  was 
an  order  of  magnitude  smaller  than  that  of  the  single-vertex  contribution.  It  was 
also  shown  that  higher -order  processes  will  dominate  the  single-vertex  diagram 
at  some  larger  value  of  n .  A  complete  report  of  pressure-induced  optical  dis¬ 
tortion  and  the  effects  of  pressure  on  the  figures  of  merit  of  candidate  materials 

was  given.  . 

^  ' 

In  the  present  report,  new  results  on  the  multiphonon  infrared  absorption 

are  reported  (SecTB),  a  theory  of  infrared  absorption  by  macroscopic  inclusions 

"  *\ 

in  crystals  is  presented  (Sec.~CJ,  experimental  curves  of  0  vs  go,  collected  from 
the  literature,  are  included  (Sec.  D),  and  a  list  of  references  to  previous  multi - 
phonon  calculations  is  given  /iirSec7  E. 

A, 

In  addition  to  this  final  technical  report,  the  following  publications  and 
conference  presentations  were  prepared  under  the  contract: 
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M.  Sparks,  "Recent  Developments  in  High-Power  Infrared  Window  Research,  ” 

4*^  ASTM  Damage  in  Laser  Materials  Symposium,  Boulder,  Colorado,  June  14-15, 
1972. 

M.  Sparks  and  M.  Cottis,  "Pres sure -Induced  Optical  Distortion  in  Infrared 
Windows, "  Conference  on  High  Power  Infrared  Laser  Window  Materials, 

Hyannis,  Massachusetts,  Oct.  30-Nov.  1,  1972  ;  proceedings  to  be  published  by 
Air  Force  Cambridge  Research  Laboratories,  Cambridge,  Massachusetts. 

M.  Sparks  and  L.  J.  Sham,  "Theory  of  Multiphonon  Infrared  Absorption, "  Con¬ 
ference  on  High  Power  Infrared  Laser  Window  Materials,  Hyannis,  Massachusetts, 
Oct.  30-Nov.  1,  1972;  proceedings  to  be  published  by  Air  Force  Cambridge  Re¬ 
search  Laboratories,  Cambridge,  Massachusetts. 

M.  Sparks  and  M.  Cottis,  "Pressure -Induced  Optical  Distortion  in  Laser  Windows, " 
to  be  published  J.  Appl.  Phys.  January  1973. 

M.  Spares  and  L.  J.  Sham,  "Exponential  Frequency  Dependence  of  Multiphonon- 
Summaaon  Infrared  Absorption,"  Solid  State  Comraun.  U,  1451  (1972). 

M.  Sparks,  "Theory  of  Infrared  Absorption  by  Macroscopic  Inclusion  in  Crystals," 
to  be  published. 

A.  Karo,  M.  Sparks,  and  L.  J.  Sham,  "Numerical  Analysis  of  Infrared  Multi¬ 
phonon  Absorption, "  in  preparation. 

Section  letters,  which  are  displayed  at  the  top  of  each  page,  are  not  dis¬ 
played  in  equation  numbers,  figure  numbers,  and  table  numbers.  Important 
results  are  denoted  by  underscored  equation  numbers. 
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B.  THEORY  OF  MULTIPHONON  INFRARED  ABSORPTION 

The  recently  observed  exponential  frequency  dependence  of  the  infrared 
absorption  coefficient  /3  in  NaCl-  and  C3  F2  -structure  crystals  is  explained  in 
terms  cf  multiphonon  absorption  processes.  A  calculation  gives  the  correct 
frequency  dependence  and  approximate  magnitude  of  0  and  predicts  the  tem¬ 
perature  dependence,  which  should  be  observable.  In  the  exponential  region, 
the  difference  processes,  in  which  some  thermally  excited  phonons  are  annihi¬ 
lated,  are  negligible  with  respect  to  the  summation  processes.  The  higher- 
order  terms  in  the  perturbation  expansion  are  negligible  for  small  n  (small 
number  of  phonons  created).  An  explanation  involving  finite  phonon  lifetimes 
is  proposed  to  explain  the  fact  that  the  alkali  halides  show  less  structure  in 
the  /3  -00  curves  than  do  the  semiconductor  crystals.  Use  of  the  central-limit 
theorem  reduces  the  mathematical  complexity  of  the  analysis  sufficiently  to 
yield  a  simple  closed- form  approximation  to  /3.  Previously  obtained  results 
that  were  given  in  the  two  preceding  reports  were  outlined  briefly  in  Sec.  A. 


I.  Introduction 


The  availability  of  high-power  infrared  lasers  has  stimulated  renewed 
interest  in  infrared  absorption.  The  intensity  I  of  infrared  radiation  propaga¬ 
ting  through  a  solid  typically  decays  according  to  the  Beer  law,  I  =  Igexp  (-0z), 

where  ,6  is  the  optical  absorption  coefficient  Materials  with  very  small  values 
-4  -1 

of  P ,  say  P  ^  10  cm  ,  are  needed  in  high -power  infrared  laser  systems 

for  such  transparent  optical  components  as  windows.  Prior  to  the  availability 

of  high-power  lasers,  there  was  considerable  interest  in  relating  the  positions 

3a 

of  the  multiphonon  peaks,  chiefly  in  the  two-phonon-summation  region  (where 
P  »  1cm  *),  to  the  phonon  dispersion  relations.  Current  interest  is  in  the 
higher -order  phonon  regions  where  P  «  1  cm  1  and  in  the  numerical  value 
of  P ,  not  simply  the  positions  of  the  multiphonon  peaks. 

2  4 

An  interesting  recent  result  ’  is  the  observation  of  optical  absorption 

( 

coefficients  that  vary  exponentially  with  frequency  j 

P  ~  exp  (-Aw)  (1.1) 

for  p  ^  10  cm"1  in  LiF,  NaCl,  KC1,  KBr,  MgF2 ,  CaF2,  BaF2 ,  SrF2 , 

Al2Og ,  and  BaTiO^ .  For  example,  P  decreases  exponentially  with  frequency 
for  over  four  orders  of  magnitude  change  in  the  value  of  P  for  NaCl  as  co/co  f 
increases  from  2.2  to  5.  8,  where  c o*  is  the  fundamental-mode  frequency  (i.  e. , 

JL 

the  frequency  of  the  transverse  optical  phonon  with  k  ^0),  See  the  points  and 
crosses  in  Fig.  1. 

The  exponential  behavior  naturally  suggests  the  form  p  ~  exp(-ftw/kgT), 
where  co  is  the  photon  frequency  and  the  other  symbols  have  their  standard 
meanings.  However,  the  room -temperature  values  of  A  in  (1. 1)  differ  by 


/ 
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•  EXPERIMENTAL  VALUES  FROM  HORRIGAN 
AND  DEUTSCH  (1972) 


•  S.  +  EXPERIMENTAL  VALUES  FROM  A.  MENTZEL, 

*  \  Z.  PHYSIK  88 ,  178  (1934). 


X. 

V 


THEORY  WITH  ONE 
'  PARAMETER  ADJUSTED 


THEORY 


\  • 
\ 


2 -phonon  3-phonon  4-phonon 

egion  region  region  / 


NaCl 


=  3. 09  x  10  rad/sec 


fi  =  L  B 


5 -phonon 


6 -phonon  l\ 
region  I 


Fig.  1.  Comparison  of  experimental  and  theoretical  results 
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factors  of  ~  2  -  4  from  the  value  of  fi  /  300  kg .  Furthermore,  the  temperature 
dependence  of  fi?  although  not  well  established  at  present,  apparently  is  less 
strong  than  exp(-fi<jO/kgT). 

Of  the  many  intrinsic  and  extrinsic  sources  of  absorption,  the  simple 
n -phonon  summation  contributions  to  fi  give  rise  to  a  frequency  dependence 
of  fi  that  is  very  nearly  exponential  over  a  large  range  of  fi ,  including  the 
experimental  range.  In  the  n-phonon  summation  process,  shown  in  Fig.  2(a), 
a  photon  is  annihilated,  a  virtual  fundamental  phonon  is  created,  the  fundamental 
phonon  is  annihilated,  and  n  phonons  are  created.  Classically,  the  electromag¬ 
netic  field  drives  the  fundamental  mode  (off  resonance,  in  general),  whose  re¬ 
laxation  frequency  is  determined  by  the  process  of  splitting  into  n  other  lattice  - 
vibrational  modes. 

By  energy  conservation,  fio!  is  equal  to  the  sum  of  the  energies  of  the  n 

final-state  phonons.  Thus,  the  n-phonon  summation  process  cannot  contribute 

to  fi  when  qj  is  greater  than  n  times  the  greatest  phonon  frequency  . 

Furthermore,  for  cc  «  noo.  ,  the  contribution  fi  of  the  n-phonon  summation 

top  n 

process  to  fi  is  small  because  energy  conservation  forces  the  wavevectors  of 

2 

the  final-state  phonons  to  be  small,  and  the  factors  k  in  the  densities  of  states 
are  small.  Thus,  fin  has  a  peak  at  a  frequency  below  nw^ .  As  n  increases, 
the  peak  shifts  to  higher  frequencies  and  decreases  in  magnitude  (since  the 
higher-order  processes  contribute  less  to  the  value  of  fi ).  This  behavior  of 
the  fin  is  illustrated  by  the  light  curves  in  Fig.  1  and  is  discussed  in  detail  in 
Sec.  IV.  The  heavy  curve,  which  is  the  sum  of  the  individual  fin ,  shows  the 
near -exponential  behavior  of  fi . 
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k 


It  will  be  shown  (Sec.  VI)  that  for  to  ^  w  at  room  temperature,  the 
contribution  to  /3  from  the  confluence  processes  (Fig.  3b  and  c)  is  negligible 
with  respect  to  that  from  the  splitting  process  (Fig.  3  a).  Assumptions  and  ap¬ 
proximations  made  in  the  analysis  are  discussed  in  Sec.  II . 

th 

In  Sec.  II ,  a  general  expression  for  the  nu  order  absorption  coefficient 

is  derived.  In  Sec.  Ill  this  expression  is  evaluated  approximately  for  n  =  2 

to  obtain  a  closed-form  approximation  for  ^  resulting  from  the  splitting  of  the 

fundamental  phonon  into  one  longitudinal  optical  and  one  longitudinal  acoustical 

mode.  In  Sec.  IV,  the  central-limit  theorem  is  used  to  reduce  the  problem  of 

evaluating  a  multiple  integral  over  phonon  coordinates  Qj  *  *  *  to  that  of 

evaluating  integrals  over  a  single  phor  on  coordinate  Q.  Further  approximations 

are  made  in  order  to  obtain  numerical  values  for  B  .  In  Sec.  V  the  central- 

n 

limit  approximation  is  discussed.  In  Sec.  VI  the  effect  of  the  confluence  pro¬ 
cesses  on  fi  is  shown  to  be  negligible.  In  Sec.  VII ,  the  effect  of  terms  other 
than  the  (n+ 1 )  -order  interaction  V^n+  ^  and  5Cn  are  discussed,  hi  Sec.  VIII, 
terms  higher  than  second  order  in  the  perturbation  expansion  are  shown  to  be 
negligible  for  small  n .  Preliminary  results  of  applying  the  Liouville  equation 
for  the  classical  distribution  function  of  an  array  cf  anharmonic  Einstein  oscil¬ 
lators  to  the  infrared  absorption  problem  are  given  in  Sec.  IX, 
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IL  Analysis  of  Absorption  Coefficient  0n 


Consider  the  an  harmonic -potential  induced  n -phonon  summation  process 
(Fig.  2  a).  It  is  not  difficult  to  show,  by  a  Green's  function  analysis  of  the  dielec¬ 
tric  constant  for  example,  that  the  value  0n  of  0  in  the  n-phonon  region  for 

£ 

materials  with  the  NaCl  or  ZnS  structure  is  [with  2r(qjj',  Cl)  in  Cowley's 

notation  written  as  T  1 
n J 


A 


n 


4ffNe*^ _ ^f  rn _ 

c  m_  n_  Cl  9  0  2  2 

*  (of- to/>  +  (Wfrn) 


(2.1) 


where  N  is  the  number  of  unit  cells,  0  is  the  volume  of  the  crystal,  e*  is  the 
Bom  effective  charge,  c  is  the  velocity  of  light,  mr  is  the  reduced  mass  of  the 
two  ion  types,  nr  is  the  index  of  refraction,  cof  is  the  fundamental  mode  fre¬ 
quency,  and  Fn  is  the  energy  relaxation  frequency  of  the  fundamental  mode. 

The  real  part  of  the  self  energy  is  not  included  in  (2. 1)  as  it  is  not  of  interest 
here;  the  experimentally  determined  phonon  frequencies  will  be  used. 

In  the  analysis  of  Tn ,  the  following  assumptions  and  approximations 
are  made: 

1)  Standard  second-order  perturbation  theory  is  used. 

2)  The  anharmonic -potential  interaction,  and  not  the  Lax-Bur stein-Bom 

7 

higher -order -moment  mechanism,  is  considered. 

3)  In  the  anharmonic  interaction,  only  the  short-range  repulsive  force 
of  nearest  neighbors  is  included. 

4)  The  central-force  approximation  is  made. 
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5)  Only  NaCl -structure  materials  are  considered.  Materials  having 
the  cubic  ZnS  structure  will  be  considered  in  a  future  report. 

6)  The  finite  lifetimes  of  the  interacting  phonons  are  neglected  except  in 
plausibility  arguments  concerning  the  sharpness  of  multiphonon  absorption 
peaks. 

7)  The  central-limit  theorem  is  used  to  reduce  the  multiple  sum  over 
phonon  coordinates  Qj  •  •  •  to  single  sums  over  the  coordinates  Q  of  one 
phonon.  In  the  present  analysis,  considerable  care  is  required  in  applying  this 
approximation,  as  discussed  in  Secs.  V  and  VI. 

8)  In  Secs.  Ill  and  IV,  approximations  to  the  density  of  phonon  states 
g(co),  the  occupation  numbers  ng,  and  certain  angle  factors  U  are  made. 

The  value  of  can  be  calculated  from  the  standard  perturbation -theory 

O 

expression  for  the  n-Boson-summation  relaxation  frequency 

T  =  i^n!  L  |vf  o  ...0  |2  A(K)6(S)n,  (2.2) 

n  Q,...Qn  »  f'Ql  n  n  n 

where  K  =k,  +  k„  +  *"+k  +  K.,k.  is  the  wavevector  of  the  i^1  phonon, 

~n  ~1  ~n 

K .  is  a  reciprocal  lattice  vector,  A  is  the  Kronecker  delta,  S  =  co  -uj~  -  •  •  •  , 

n  «^l 

Q.  is  defined  as  the  set  of  numbers  k .  and  b. ,  where  b.  is  the  branch  of  the  dis- 
i  '“il  l 

persion  relation,  =  (nj  +  1X^  +  1)  ••  •  (nn+ 1  )/(nw  +1 ),  nj  =  n(Q. ) 

=  [exp(o)g.  /wT)  -l]"1,  n^  =  [exp(w/toT)  -  l]"1,  coT  =  kgT/f^  3.928X  1013 
rad/sec  at  300  K,  and  the  V's  are  the  coefficients  in  the  scattering  Hamiltonian 


K  =  £  V.  n  n  A(K)  a.an  •  •  •  an 

Ql*',Qn  f’Ql’*,Qn  ~  f  Ql  Qn 


+  cc  . 


(2.3) 
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In  (2. 3),  f  denotes  the  fundamental  mode,  cc  denotes  hermitian  conjugate,  and 
the  a"^"  and  a  are  phonon  creation  and  annihilation  operators  normalized  to  unit 
commutators. 

For  the  Hamiltonian,  up  to  nul  order  phonon  processes,  take  the  central- 
force  potential 


K  = 


1 

fc+ry: 


,(nfl) 


E 

A= 1 


N-l 

E 

m=l 


^m  ®£m 


n+1  n 

)  +  E 

p=2 


/(p) 


Am 


(2.4) 


I 

i 

i 


where  R^m  is  the  length  of  the  vector  from  ion  m  to  ion  A  ,  d^n+^V(R^  )/ 

dRJjrn^n+^  evaluated  at  equilibrium  (R^m  =  6^  ),  A  labels  the  cell  (lighter  ions, 
for  example)  and  m  labels  the  neighboring  ions.  No  factor  of  1  /2  to  avoid  double 
counting  is  needed  for  this  labeling.  The  vector  amplitude  R^m ,  whose  equilibrium 
value  is  6^m ,  can  be  expanded  in  terms  of  displacements  u^  and  um  as 


-  R  _  +  u. 
'-mO  l 


u 

~m 


^  ~im  +  ~£m 


9 


where  JR^  and  RmQ  are  the  equilibrium  values  of  the  positions  R^  and  Rm  of 
the  and  mth  ions  and 


6 .  =6  =  R._  -  R  , 

~jsm  «m  ~j20  ~m0 


Si  =  Ji%0>  * 


i  =  A  ,  m. 


~2m  ~m 


The  amplitude  R^m  can  be  written  as 

Rim  ■  NSm  +  W  •  <£m  +  2-m)] 


1/2 


=  6  +6  •  u .  + 

m  m  m 


•  »  •  • 


(2.5) 
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R„  -6  =6  •  u .  + 

£m  m  m  ~im 


Dropping  the  p  =  1  •  •  •  n  terms  in  (2.4)  and  the  higher -order  terms  in  (2.6)  and 
substituting  (2. 6)  into  (2. 4)  gives 


1  vOrt-D 


K„  *  (Kmrv  £“«•*■«> 


When  all  terms  from  (2.4)  and  (2.6)  are  retained,  the  angle  factor  6m  •  u^m  in 
(2.7)  [and  an  in  (2. 16)]  is  changed,  as  discussed  in  Sec.  VII. 

9 

The  usual  expansion  of  the  displacement  is 


/  *  \1/2 

ttT<BiO>  *  *  (  ?FTmTaQ  )  HtQ  Aq 


ik  •  R.n 

**  ~i0 


where 


AQ  "  \b  ~  \b  +  a-kb 

and  r  =  1  or  2  denotes  the  type  of  ion  (Na  or  Cl,  for  example).  The  polarization 
vectors  utq  obey  the  orthonormalization  relations 


h  u  , .*  •  u„, ,  ,  =  6,  u, 
~7\kb  ~7kb  bb 
T 


(2.9a) 


(2. 9b) 


where  x^  is  a  unit  vector  along  a  unit- cube  axis.  These  are  simply  the  ortho¬ 
normality  relations  for  the  rows  and  columns  of  the  six-by-six  diagonalizing 
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matrix  (for  the  eigenfrequency  problem)  whose  columns  are  the  six  vectors 
(u>o*  £<q)*  where  T  =  >  and  T  =<  denote  the  greater  and  lesser  masses, 
respectively. 

In  Ref.  9,  (2. 8)  was  written  [with  x  -  T ,  j  =  b,  (b ,  k  )  =  Q,  and 
u(ix)  =  ut(RX0)]  as 

1/2  i 

=  q  (  2NmTOJQ  )  £<TUb>  AQe 


ik  •  R  cr. 


In  this  expression,  the  equilibrium  ion  position  appears  in  the  exponent,  in 
contrast  to  (2.8)  where  R.n  =  Rtfft  +  B.  (with  B.  as  one  of  the  two  basis  vectors 
specifying  the  ion  positions  in  a  cell)  appears  in  the  exponent.  The  relation  be¬ 
tween  these  two  expressions  is 

utq  =  e(rjkb)  ,  for  small  mass, 

ik  *  B 

JiTQ  =  ~(Tl-b)e  ^  *  for  large  mass. 


where  it  is  assumed  that  the  small  masses  are  at  the  lattice  sites  ( B  j  -  0 )  and 
the  large  masses  have  Jg2  H  B  =  a^  (x  +  y  +  z)  for  the  NaCl  structure.  Here 
denotes  a  unit  vector  and  ann  is  the  near  neighbor  equilibrium  spacing. 

For  the  fundamental  mode 

1/2 

(2. 10a) 


T  /  m  \  ' 

««  =  <-*>  (57)  sf 


where  mr  is  the  reduced  mass  (m^.'1  =  m^1  +  m^1 ).  In  the  limit  kL  «  1 
where  L  is  of  the  order  of  the  unit  cell  length. 
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T  /">r\  - 

~Tkb  ^  ^  (  mT  )  ^Tkb  ’ 

,  1/2 
/  mT  \  - 

~Tkb  "  \  m<  +  m>  /  UTkb  ' 


(optical) 


(acoustical)  . 


(2. 10b) 


(2. 10c) 


Notice  that  (2. 10a)  and  (2. 8)  with  /=  0  and  all  other  Aq  =  0  (only 
fundamental  phonons  present)  give  the  well  known  relations  ^ 


Br(Si>  -  £t'  (£i  +  ~m }  ~  1//7^ 


u_(R.  )/u_,  (R.  +6  )  ~  -m  ./m_ 

Tx~i"  r  ~i  ~m  T  T 


for  T  £  T'. 


Substituting  (2.  8)-(2. 10a)  into  (2.7)  and  using 


L  exp[i(ki+...+kn)  •  Rxq]  =  NA(Kr' 
X — 1 


Him  5  -  am  5  a<<Si>  -  2J><5io +  6m> 


/  *  \V2u 

[  TN  "m<ooQ  )  U 


m(Q>  AQe 


where 


U  (Q)  =  u -  /  -  u ^(~ 

~mv  ~<Q  sj  m>  ~ >Q 


ik  •  6 


(2. 11) 
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gives  (2.3)  with 


l.v  iGH-1) 


(■jh)  7  v(n+1Wl)  i  i  i 

fQl,,,Qn'  i(n-l)  ,  .1/2  m  n/2  1/2 


X  £  <6m‘VUm^l>-Um<Qn> 


(2. 12) 


where  Um(Q)  =  6m  *  In  general. 


6  •  U  ,(Q)  5  U  ,(Q)  , 

m  ~m  m  ' 


Um«»  -  Um(Q) 


(2. 12a) 


The  factor  (n+1)  in  (2. 12)  arises  from  the  n+1  ways  of  choosing  a^  from  the 
Aq.  in  (2.7)  and  (2.  8).  Note  that  for  the  fundamental  mode,  (2. 11)  and  (2. 10a) 


.  m  .  1/2 

Bm™- (.7)  ai  • 


(2. 12b) 


The  following  model  potential  *  is  used:  V  =  E  exp(-R£m/pa  )  for  near¬ 
est  neighbors  and  zero  otherwise.  Here  a^  is  the  equilibrium  near -neighbor 
distance.  The  Coulomb  potential  is  added  to  this  repulsive  potential  in  order  to 
determine  the  values  of  p  and  of  E  =  3exp( l/pJBa^  p  f(l  -  2p),  where  B  is 
the  bulk  modulus,  but  is  not  included  in  the  anharmonic  potential.  This  is  a  rea¬ 
sonable  approximation  for  2  ^  n  ^  20,  since  the  higher -order  derivatives  of 
the  rapidly  varying  repulsive  potential  are  much  greater  than  those  of  the  slowly 
varying  Coulomb  potential.  This  is  easily  seen  for  a  repulsive  potential  ~1/R^ 

Hi 

and  a  Coulomb  potential  ~1/Rjj »  for  which  the  n  derivatives  contain  factors  of 
magnitude  (2)(3) *  -  *  (n+1)  and  (10)(ll)",(n+9),  respectively. 
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The  value  of  V'  1  for  this  model  potential  is 


3Ba  zp  , 

nn  1 

d-2P)  nn 

nn  H 


(2. 13) 


From  (2. 12b) 


1/2  „ 

U  (f)  =  (nv,/m  )  6  •  u, 

nr  '  < '  r '  m  f 


(2. 14) 


Substituting  (2. 13)  and  (2. 14)  into  (2. 12)  gives 


3  P  ® ann  | 

•Q„=  ±JT7W)  °n  n,N/(n-l)  ^-—TJTT 

^1  ^n 


\  1/2  i 

Jl _  \  f  ft 

"%/  [2p2  m<  _ 


(2. 15) 


where 


mE,  <V  °m«l> 

m=l 


U  (Q  ) 
m  n 


(2. 16) 


Substituting  (2. 15)  into  (2.2)  gives 


r  =  y  Z  , 

n  n  n 


(2. 17) 
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where 


x  [l-exp(-ti<J/kBT)] 


and 


Z  = 


_L 

N 


Z 

Q, 


)+  l 


a>. 


1 

TT 


n(a,Qn  )+1 


\2 

I  N  A(£n>  6  ("  -  wQi . 0>Q  )  . 


Substituting  (2. 17)  into  (2. 1)  gives,  for  O)^  »  0)  ^  +  T  ^ 


0 


S'  U  Ku 

2  p  mx  n, .  + 


03 


1 

0)4n! 


(03  D  )  Z 
mx  e'  n 


(2. 18) 


(2. 19) 


where 


(noj  +  1)  =  1  -exp  (-fiOJ/kgT) 


— 2 — 5 - 

2p  ann2  m. 
nn  < 


uj 

mx 
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3* 


and  w  ,  which  cancels  out  of  (2. 19),  was  introduced  for  later  convenience. 
From  (2. 18)  and  (2. 16), 


Sn  = 


6 

E 

m=l 


6 

E 

m'=l 


(6  *uf  )(6 


m 


m 


uf) 


1 

TJ 


E 

Ql 


n(ajQi)+i 


tc 


Qi 


Um«VUm*«V 


1 

TI 


E 

Q 


n<“Qn)+1 


n 


cc 


Qn 


U  (Q  )  U  , 
m  n  m 


«„> 


x  NA(Kn)6(o,-idQi . wQn)  • 


(2.20) 


A 

It  can  be  shown  that  /3n  is  independent  of  the  direction  of  u^  for  the  NaCl 
structure.  The  proof  involves  using 


U  (b.,-k. )  =  U  (b. ,  k. ) 

m'  i’  ~i '  m v  l  ~i 


(2.21) 


and 


U  (Q.)  =  -U  (Q. )  , 

m'n7  -m  l 


(2.22) 


where  U_m(Q.)  is  obtained  from  (2. 11)  and  (2. 12a)  by  replacing  by  -6m. 
These  results  (2.21)  and  (2.22),  which  follow  directly  from  (2. 11)  and  (2. 12a), 
are  used  in  (2. 20),  the  dummy  indices  in  the  sums  over  the  components  of  the 
wavevectors  are  relabeled,  and  the  invariance  of  Wq.  to  the  symmetry  operations 
of  the  cubic  crystal  are  used  to  complete  the  proof. 

Since  the  value  of  0n  is  independent  of  direction  of  u^  ,  we  choose 

uf  =  Xj  ,  (2.23) 


/X 

where  Xj  is  a  unit  vector  along  a  unit-cube  axis.  Then,  only  the  two  terms  with 

6  =  ±  a  x.  in  (2. 16)  are  nonvanishing: 

~m  nn  1  ° 


-  ux(Qi)-ux(Qn)-ux(QJ)-ux(Qn)  . 


(2.24) 
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III.  Closed-Form  Approximation  to  f>„ 


For  n  =  2  ,  (2. 24)  gives 


a2  = 


<QjX 


<Q2x 


-  I^<  u  eiannklx\ 

v®>  u>QjX  e  I 


r™< 

j  m> u 


>Q2x 


■iannklx 


(“<V  '  U>Q]^  lam,,t1^ 

("<Q2x  -  u>Q2x  j 


(3.1) 


where  subscript  x  denotes  the  Xj  component,  and  wavevector  conservation. 
k2x  =  “  klx’  has  been  used*  Combining  terms  in  (3. 1)  and  using 

~T  -  kb  ~Tkb 


gives 


a2  =  2i 


/ 

Vm>  (u<1cIbIxu>kIb2x-“>k1b1x“<k1b2x)  8i"llx‘ 


From  (2. 18)  with  n  =  2 


£  _  JL_  j,  n(wklb2)+i 

2  ‘ F  iibib2  Vi  ~V7~ 


“2  I  6(“-‘Vr“k,b2)  • 
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The  value  of  a2  tends  to  be  small  in  a  number  of  cases.  First,  for  both 

phonons  (Qj  and  Q2)  on  the  same  branch,  i.e. ,  bj  =  ,  (3.2)  gives  the  well 

7  11 

known  selection  rule  *  =  G;  i.  e. , 

P2  =  0  ,  for  b1  =  b2  .  (3.4) 

I  2 

For  kja^  «  1,  |  sin  k^a^  -  (^ixann)  <<c  **  and  the  contribution 
to  ^  is  small. 

For  the  one-dimensional  diatomic  chain  (Ref.  10,  p.  176  ,  for  example) 
and  for  the  three-dimensional  NaCl  structure  for  k  along  [  111  ] ,  the  optical 
and  acoustical  modes  at  the  zone  boundary  correspond  to  standing  waves  with 
only  the  lighter  and  heavier  masses,  respectively,  vibrating: 


Et<Q  =  0  ’  ~r>0  =  Gt>0  ’  acoustical,  (3.5) 

x  x  x  zone  edge 


A 

~T<Q  =  UT<Q  ’ 


ut>Q  =  0  ' 


optical,  (3. 6) 

zone  edge  . 


If  both  modes  in  (3. 2)  obey  either  (3. 5)  or  (3. 6),  then  a2  =  0 .  Even  though 

(3. 5)  and  (3. 6)  are  not  satisfied  exactly  for  other  points  in  the  zone,  this  result 

suggests  that  the  quasi- selection  rule  that  ftp  is  small  unless  one  output  phonon 

is  optical  and  the  other  acoustical  may  hold  for  NaCl -structure  crystals.  The 

validity  of  this  proposed  quasi- selection  rule,  which  has  several  important  con- 
3 

sequences,  will  be  examined  in  detail  in  a  later  report. 

The  sum  in  (3. 3)  can  be  approximated  as  follows.  Consider  the  term 
bj  =  LO  and  b2  =  LA  and  the  term  bj  =  LA  and  b2  =  LO ,  where  LO  and  LA 
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denote  longitudinal  optical  and  longitudinal  acoustical,  respectively.  Since  (3. 3) 
is  unchanged  if  bj  and  b2  are  interchanged,  these  two  terms  give  the  same 
contribution  to  £2  • 

i  .2 

A  factor  of  1/3  will  be  included  in  the  value  of  |  a2  i  to  account  for 

i  i2 

several  regions  in  k  space  in  which  \ol^  |  is  small.  These  regions  are  as 
follows:  The  factor  sink^a^  in  (3.2)  makes  the  contribution  to  S2  from  the 

i  |2 

volume  in  k  space  near  the  plane  k^x  =  0  small.  Similarly,  |a2  I  is  small 

for  all  directions  when  kann  <<  1 .  At  the  zone  edges  along  [  100]  axes, 

ka  =  ir,  and  the  factor  sink,  a  in  (3.2)  is  equal  to  zero, 
nn  lx  nn 

,  ,2 

In  the  rest  of  k  space,  the  value  of  |  I  will  be  approximated  by  the 

value 


2  _  4 

111,  ZB 


9 


(3.7) 


at  the  zone  boundary  along  a  [  111  ]  axis.  [  Note  that  at  this  symmetry  point, 
which  is  commonly  denoted  L,  (3.5)  and  (3. 6)  give 


u>  LO  "  U<LA  =  0  ’  u<LO  '  U>LA  '  1  * 


(3.8) 


Since  the  x  component  of  a  unit  vector  along  the  [  111 )  axis  is  l/ZT1,  the  value 

of  the  square  of  the  u  factor  in  parentheses  in  (3.2)  is  1/9.]  With  the  factor-of-3 
2 

reduction  in  |  or2  |  discussed  in  the  preceding  paragraph,  and  a  factor  of  two 
included  for  the  two  choices  of  bj  and  b2  discussed  above. 


m. 


ITU 


f. 

•1 

4 

i 


8 

=  Yf 


(3.9) 


aec.  a 


The  density  of  longitudinal  acoustical  states  is  rather  sharply  peaked  atJ 


=  2.7  x  10  6  rad/sec. 


(3. 10) 


(See  Fig.  4. )  Thus,  to^  will  be  set  equal  to  the  constant  to  ..  .  With  these 

~  1 

approximations,  (3. 3)  and  (3. 9)  give 


o  n(to  )+l  n(to  )+l 

'2LL  =  Tl  S7  -s^~  v  ks  —sr  6  [V<“-“la>]  <3-"> 


for  the  contribution  to  from  t^ie  LO  and  LA  modes.  Since  the  summand 
now  is  a  function  of  to  ^  only,  the  sum  on  k  j  can  be  replaced  by  an  integral 


over  to . 


S  f(tOj)  s  ]  dtOj  g^to^fftOj)  , 

k.  Jn 


where  the  density  of  optical  states  gQ(to)  is  normalized  to  unity 


dto  gQ(to)  =  1 


Evaluating  the  trivial  integral  in  (3. 11)  gives 


„  8  m<  . .  "<“LA>+1  "'“•“’lA>+1 

E2LL  “  27  m^g0(“'“LA*  £T7  <0  -<0.  , 

>  LA  LA 


(3. 12) 


(3. 13) 


A  rough  sketch  of  ^ ^  from  (3. 13)  and  (2. 17)  is  shown  in  Fig.  5  as  a  dashed 
curve.  It  will  be  argued  below  that  the  effect  of  replacing  the  delta  function 


27 


X-^Ski.u t 
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for  the  value  of  ^LL  at  t^ie  Pea^  at  u  ~  6.  3  x  10  rad/sec.  In  units  of 


1013  rad/sec 


WLA  2,7  ' 


to  =  3. 6  , 
gP 


“>LO  =  '"9  • 


to. +  to  =  6.  3  , 
LO  gp 


toT  =  3.93  , 


.12 


at  T  =  300  K  . 

Substituting  (TOI .  )  =  1.0  x  10A*  rad/sec  from  (3. 15)  into  (2. 1)  gives 

Z  LL  my 

(3.  16) 


(/39)  =  19cm 

“  LLmx 


mx 

1 


A  comparison  of  this  result  with  the  numerical  results  currently  underway  and 
a  comparison  of  the  numerical  results  with  experimental  results  will  be  given 
in  a  future  report. 

A  second  approximation  to  the  value  of  £2LL  can  be  obtailied  approxi¬ 
mating  the  Brillouin  zone  by  a  sphere  of  radius  k^  and  selecting  linear  disper¬ 
sion  relations  for  the  LO  and  LA  branches  to  match  the  densities  of  states. 

For  the  LA  branch,  take 


"kLA  =  vA(,t-W+"A- 

See  Fig.  6a .  Then  using 

k 


VA  5  ^A’^^mx 


(3. 17) 


k 

(V 


3  n  r  ™ 
mx  JQ 


dkk*  f(k) 


(3. 18) 
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gLA<“>  *  l  6  <“iLA  -  “> 


gives 


*la(">  = 


dkk  6  (co.  -  co) 


(3. 19) 


Using  (3. 17),  the  delta  function  can  be  written  as 


6<cuVt  .  -  co)  =  6[vA(k  -  kfflX)+  coA  -co]=  6(k  -  kA )  , 


where 


kA  "  kmx  -  "7 


coA  -  co  co  -  coA 


<V*a 


Using  these  last  two  results  and  evaluating  the  integral  in  (3. 19)  gives 


3N«o-wa)2 

giA(w)  =  - “ ~~T  » 


for  coA  <  co  <  coA 


(3.20) 


=  0  , 


otherwise  . 


For  the  LO  branch,  take 


“V-LO  =  -V0(k  ‘  kmx)+  “d  ’  v0  =  <<V“o)/km*  • 


(3.21) 


See  Fig.  6a.  The  same  type  of  analysis  as  above  gives 


3N(co  -cof 

gLO(w)  =  - - T 

<"o  -  “o> 


for  coQ  <  co  <  coQ 


(3. 22) 


=  0  , 


otherwise. 
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Hie  densities  of  states  in  (3.20)  and  (3.22)  are  sketched  in  Fig.  6b.  Formally 
inserting  a  sine  factor  [  see  (3.2)]  in  (3. 9)  gives 


a  |2_  8  mNa  .  2, 

2  !  ~  2 T  ~  sin  k’  a 


Cl 


1a  • 

nn 


With  (3. 3),  this  gives 


2LL  "  27 


8  mNa  3  f  1 

27  mCl  k  3  ] 

not  Jg 


”“dkk2  "(“k^)+I 


03, 


CLA 


n(“kLO)+i 


03*.  T  ~  ^  (k  -  k.  )  sin^  k,  a  , 

KLO  L'  1  nn  ' 


(3.23) 


where  03.  .  and 


kLA  and  ^LO  3re  817611  by  (3‘  17>  and  (3.21),  respectively,  and 


k  -  “0^-“ 


W0+WA-(O0  -03a  ^ 
Evaluating  the  integral  in  (3.23)  gives 


n  (“kL  LA )  +  1 


r  -  A  mNa  H  + 

2LL  ~  T  THT"  “= - = - g  - - - 

Cl  (%  +  <VV“a> 


"H  LO)+1 

KLLU  .  2 . 

77 - —  sin  k.  a  , 

\LO  L  nn 


for  U30+  03a  <  03  <  U3q  +  U3a 


=  0 

Here  03.  T  .  and 


otherwise  . 


(3.24) 


kLLA  and  wkLLO  are  the  valuaa  of  ^  ^  and  03^  from  (3. 


(3.21)  with  k  =  kL. 


17)  and 
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The  maximum  value  of  (3. 24),  which  occurs  at  to  =  tOg  +  u.'A  (i,  e. , 


kL  =  km*>-is 


8  m< 


n(wA)  +  l 


>  +  *WA> 


n(w0)+  1 


(3.25) 


where  sin  ann  was  set  equal  to  one  formally.  In  units  of  10  rad/sec. 


WLO  4,9  ' 


U)n  =  oj  =  3. 5  , 
0  gP 


wA  =  2.0  , 


u>A  =  2.8  , 


%  +  WA  =  6'9  * 


coQ  +  ioA  =  6.3  . 


Thus 


W0  +  WA  "  W0  "  =  x  *0  rad/ sec  =  toLO/8.2 


and  (3. 25)  can  be  written  as 


7#  3  rn<  n(wA)+l  n(Wg)+l 


mx  WL0  m> 


(3.26) 


JCC.  L 


In  comparing  (3.26)  with  the  previous  result  (3. 14),  the  unmodified  curve 

(dashed)  in  Fig.  5  should  be  used.  Thus,  (3. 14)  should  be  multiplied  by  a  factor 

of  1. 9,  which  is  the  ratio  of  the  peaks  of  the  two  curves  in  Fig.  5.  Then,  with 

13  13 

(Cq  =  u  and  cu.  =  2. 8  x  10  rad/sec  and  2. 7x  10  rad/sec,  the 

gP 

value  of  in  (3.26)  is  approximately  a  factor  of  1.5  times  larger  than  the 
corresponding  value  in  (3. 14),  which  is  good  agreement. 

In  the  second  model  [spherical  Brillouin  zone  and  (t)  =  f(|jc  | )] ,  setting 
a>A  =  u>A  =  it) ^  corresponds  to  using  the  delta  function  approximation  to  the 
density  of  LA  states  gjj^(u)).  For  this  second  model,  it  will  now  be  shown 
that  the  effect  of  replacing  the  delta-function  approximation  to  g^A(co)  by  the 
quadratic  density  (3. 20)  (with  it) A  ^  a>A)  is  to  reduce  the  region  over  which 
T2ll  is  nonzero  and  increase  the  maximum  value  of 

Setting  5?a  =  <CA  =  <i.e.  -  «A  -  *  "A=  “0  -“>0  =  “>LO/3-5) 


in  (3.25)  gives 


mx 


3.1 

it) 


m 


<  n(WLA)+1  n^S£}+l 


LO 


m. 


it) 


LA 


it) 


(3.27) 


gP 


as  the  maximum  value  of  (3.24)  in  the  delta -function  approximation.  The  delta- 
function  approximation  to  (3. 24)  is  sketched  as  the  dashed  curve  in  Fig.  7.  The 
corresponding  curve  for  the  density  gLA(w)  in  (3.20)  is  shown  as  the  solid  curve. 


This  result  that  the  delta -function  approximation  to  ^LL  is  nonzero  over  a 
larger  region  is  surprising  in  view  of  the  result  in  Sec.  IV  that  ^LL  is  a  convolu' 
tion  of  the  two  densities  gj^  and  g^Q  (weighted  appropriately).  Recall  that  the 


< 
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Fig.  7.  Relaxation  frequency  T2LL  from  the  spherical-Brillouin-zone 
approximation  (dashed  line).  Making  the  delta  function  approximation  to 
g^to)  gives  the  solid  curve. 
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length  of  the  nonzero  region  of  convolved  functions  is  the  sum  of  the  lengths 
of  the  nonzero  regions  of  the  individual  functions.  The  difference  in  these  two 
results  arises  from  the  fact  that  wavevector  conservation  was  neglected  in 
deriving  the  convolution  result. 


This  difference  can  be  understood  in  terms  of  the  schematic  illustrations 
in  Fig.  8.  The  dispersion  curves  for  the  LO  and  LA  modes  are  shown  (for  a 
spherical  Brillouin  zone)  at  the  bottom  of  each  section  (a)  •  •  •  (d)  of  the  figure. 

The  bars  at  the  top  of  each  section  indicate  the  nonzero  regions  of  (i.  e. , 
the  region  of  uj  for  which  cc  =  to  ^  +  to  for  some  pair  of  modes.  The  bars 
to  the  left  are  for  the  case  of  =  k^Q ,  and  those  to  the  right  are  for  no  wave- 
vector  conservation.  In  going  from  Fig.  8c  [  for  which  g^(to)=  6  (to  -  to^)] 
to  Fig.  8b  [for  which  gjj^60)  is  given  by  (3.20)],  the  left-hand  bar  decreases 
in  length  while  the  right-hand  bar  increases. 
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Fig.  8.  Schematic  illustration  of  regions  for  which  the  two-phonon- 
summation  process  with  to  =  to  ^  +  to^  can  conserve  energy.  The  bars 
at  the  top  left  of  each  figure  show  this  region  for  the  case  in  which  k  is 
conserved,  and  those  to  the  right  are  for  no  k  conservation.  The  model 
dispersion  curves  for  the  four  cases  of  interest  are  shown  in  the  bottom 
portions  of  the  figures. 
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IV.  Closed- Form  Approximation  to  0n 

The  value  of  in  (2. 18)  for  general  n  will  be  approximated  by  using  the 
central-limit  theorem  as  follows:  First,  wavevector  conservation  will  be  neg¬ 
lected  by  formally  setting 

N  A(K  )  =  1  .  (4.1) 

This  is  expected  to  be  a  good  approximation  for  large  n  .  Substituting  (2.23) 
into  (2.20)  and  using  (4. 1),  (2.21),  and  (2.22)  gives 

rri-1 

E  =  E  .  +  (-1)  E  ,  (4.2) 

n  n+  '  n- 


where 


1  Ql  1 

E  =  2  ~  £  — —  0  (Q.) 
n±  TT  Jr  u:n  w±wr 
Ql  Qi 


nrv  + 1 

*  %  ^  w 


X  6(ChJ-«q . OJq^)  , 


(4.3) 


with 


0±(Q.)  h  [ReUx(Q.)]2±  [JmUx(Qi)]2 


By  inserting  the  trivial  relations 


dfif(fl)  6(0-  WQ  )  =  f(wQ  ) 


(4.4) 


into  (4. 3),  En+  can  be  written  as 


n 


n± 


2<a+>  f  ( 

-  — ir  )  dV±(0i>- )  dV±«V 

mx  -®  -oo 


n 

x  6  (to  -  E  0. )  , 

i=l  1 


(4.5) 
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where 


ff+(fl)  = 


mx  z  n(Q)+l  0+(Q) 


N  <a+>  q  Wq 


x  6  (0  -  wQ) 


and 


<“±>‘  \  d0  4  = 


(4.6) 


6  (fi  -  o)q) 


U) 


nfaJn  )+ 1 


-2HL  2  - ^ -  0  (Q)  . 

N  q  “q 


(4.7) 


Note  that  o+  (  fi )  are  normalized  to  unity: 


r 

l  dfi  CT+ 


(0)  =  1  . 


The  factor  of  ,  which  cancels  out  when  (4. 6)  is  substituted  into  (4. 5),  was 


introduced  for  later  convenience. 
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For  large  n ,  the  central -limit  theorem  reduces  (4. 5)  to 


2<ttt>D  -<“-"“±“n«>2/2"“o.2“rnx2 


n± 


flvn  a  iomv 
a+  mx 


n+1 


(4.8) 


where 


A® 

w  =  1  dfi  a  (fi)  fl  s  a.  a) 

±1  ±  ±  mx 

V-® 


(4.  9) 
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and 


03, 


■[ 


d0  o±  ( O)  n2  =  a,2  ujf 


(4. 10) 


and 


03 


2  _  2  —2  _ 


a+ 


ui  r  -  03  j.  2  tt  03 
±  ±  <t+  mx 


(4. 11) 


With  a+(  O)  given  by  (4.6),  the  results  (4.9)  and  (4. 10)  become 


ff±*  k T<k)  l  ("Q+1>^(Q) 


and 


a+2  ~ 


N(a±>umx  Q  <nQ+1)UQ  *±(Q) 


Substituting  (4. 3)  and  (4.  8)  into  (2. 19)  gives 


.n+1 


(4. 12) 


(4. 13) 


where  the  approximate  equality  is  discussed  in  the  following  paragraph  and 
.  DPK  /“ t-nu.  1  l” 

^n±  «V  ( n  + 1 )  '  “  '  /?•!  (  *  e) 

±  ff  +  '  03 


x  exp 


—  2  9  9” 

-  (oj  -  na,  co  )  /2na„  4oi_  ^ 

'  ±  mx'  '  a+  mx 


(4. 14) 


Using  the  central-limit  theorem  has  reduced  the  multiple  sum  over  Qj  •  *  [in 
(2.20)]  to  sums  over  the  single  phonon  coordinate  Q  [in  (4.7),  (4. 12),  and  (4. 13)]. 

By  calculating  the  values  of  the  real  and  imaginary  parts  of  Ux(Q.)  for  several 

14 

points  in  the  Brillouin  zone  using  the  eigenvector  data  of  Karo  and  Hardy ,  it  is 
simple  to  show  that  neither  of  the  two  parts  is  negligible  in  general.  Thus,  from 
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the  definition  of  0  +  (Q . )  under  (4.3),  0+  =  0_/s,  where  s  ^  2  .  From  (4. 3), 
-  En_  /  sn ;  thus  £n_  is  negligible  for  large  n  since  0n+  «  Lr+  . 


The  values  of  the  a's  will  now  be  estimated.  First  consider  Ux(Q) .  For 
kL  «  1 ,  the  approximations  (2. 10)  are  valid.  For  modes  near  the  zone  edges, 
the  one-dimensional  results  (3.5)  and  (3. 6)  will  be  used. 

Substituting  these  results  into  (2. 12)  gives 


A  A 


Um<Q>!  ^<6m*u<0>  • 


optical, 
zone  edge 


2  m<  -  -  2 

V®  *  -S-  <V">Q>  • 


acoustical, 
zone  edge 


U  (Q)  I  as  ^  (8  .  u  )2 

mr  '  m  tQ'  * 


m 


optical, 
kL  «  1 


A  A 


A  A 


2  2 


Um«>|  ~  (6m‘  utQ)  <V  V  <kL>  . 


acoustical, 
kL  «  1  . 


Since  kL  «  1  in  the  last  expression,  this  value  of  Um(Q.) 


is  negligible. 


Thus,  as  a  rough  approximation,  Um(Qj)  will  be  set  to  zero  for  oj  <  f 
where  f  <  1 .  In  the  nonzero  region,  the  remaining  three  terms  are  averaged 
to  give 

I  I2  ?/  \  „  2 

I  Ux<e>|  H-t(1+  S^)«V“tQ>  ’  <4- 15> 


where  6  is  the  unit  step  function  and 


*(' 


ti1  +  sn 


m<  \ _  2  ( 
m>)=  3  [ 


1  + 
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The  value  of  the  scalar  product  in  (4. 15)  is  approximated  by  1  //  3 

^  A 

which  is  its  value  for  6m  =  Xj  and  u<j  a  unit  vector  along  [111].  The 
value  of  the  square  root  of  the  average  value  of  cos^  0  also  is  equal  to  l//"^ : 


■f-jj-  IdOcos^S  =  -j-  \  d0sin0cos^0 


du  u  =  1/3  . 


This  reduces  (4. 15)  to 


|u  (0)  as  -l  (i+  —  )  6  (fa)  -  U)n)  . 

!  x  '  y  y  m>  /  '  mx  Q 


(4. 16) 


The  sum  in  (4. 7)  now  can  be  written  as 


F^F(V=  6 


dWQ  S<"q)  F<WQ> 


(4. 17) 


where  the  total  density  of  states  g(  u>)  is  normalized  to  unity 


doo  g(co)  =  1 


in  order  to  make  the  total  number  of  states  correct: 


L  (1)  =  6  N 

Q 


dco  g(co)  = 


Sec.  B 


Substituting  (4. 16)  into  (4.7)  and  using  (4. 17)  gives 

<  o.  >  =  -s-  (  i+  — -I  <  )  u 

where 


v  4  (.*£)<¥•) 


mx  ’ 


AGD 

\  dojg(«)  . 

^fco 

mx 


Similarly  (4. 12)  and  (4. 13)  give 


a, 


< (n+ 1 ) 9) 


and 


_  ( (n+l)ui6> 


+  /2±ie\u  2  ’ 

\  to  /  mx 


where 


s 


((n+l)9>  =  \  dw  g(w)  [n(co)+l] 

hui 


mx 


and 


•S' 

Vf/.t 


((n+l)w0)=  \  das  g(io)  [n(u))+l]  oo 

'too 


mx 


(4. 18) 


(4. 19) 


(4.20) 


(4.21) 
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The  density  of  states,  shown  as  the  solid  curve  in  Fig.  9,  will  be 
approximated  by  the  Debye  density 


g(«>  =  * 


(4.22) 


CO 


mx 


as  illustrated  by  the  dashed  curve  in  Fig.  9.  The  factor  n(to)+l  will  be  approxi¬ 
mated  by  its  high -temperature  limit 


n(co)  +  l  = 


1 


to. 


1 


"  -to/to™  ^  to  +  T 
1  -e  '  T 


(4.23) 


In  the  worst  case  of  to/to^.  —  1  for  NaCl  at  300 K,  n(to)+l  =  1.582  and 
(co^/co)  +  =  1.500;  thus,  the  approximation  ?r  (4.23)  is  quite  good. 


Substituting  (4.21)  and  (4.22)  into  (4. 18)-(4.21)  and  evaluating  the  simple 
integrals  gives 


(4.24) 


(Cr 


= 


2  T  1  3 

2[(1'f  >  «£  +  T(1-f  >3 


- SJ 

4(1  -f)  w 


T  ..  .2 


mx 


+  (1-0 


(4.25) 


7- 


(1-f) 


(0 


to 


-  +  4  d-f2) 


mx 


(4.26) 
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and 


6l' 


mx 


3(1  "O  ^-+  4d-f2) 

mx 


(4.27) 


Substituting  these  results  (4.24)-(4.27)  into  (4. 14)  and  assuming  that  the  con¬ 
tribution  from  the  j3n_  term  is  negligible  as  discussed  under  (4. 14)  gives 


*n  - 


D-K  /  a, to  \  . 

_ p _  I  1  mx  j  1 

a,4a  (n  +1)  '  “  '  StT  nl 

1  01  Ui 


\±  (i+m-±)  At!*,  e\  D  l"  exp 

/  -  x2 
(oe-nof.oe  ) 

'  1  mx' 

3  \  m>  /  oe  ^ mx  ?  e 

~  2  T~ 

2na„  oe 

L  J 

a  mx 

(4.28) 


where  and  K  are  defined  under  (2. 19). 


1  3  1 

Consider  two  values  of  f:  f  =  y  and  f  =  y  (i.e. ,  f  =  0.79).  The 
significance  of  considering  more  than  one  value  of  f  will  be  discussed  in  Sec.  V. 
The  value  of  f  =  y  corresponds  to  assuming  that  kL  «  1  for  1/8  of  the  modes 
(or  1/4  of  the  acoustical  modes)  according  to  the  discussion  of  (4. 15)  above. 

Note  that  with  g(uj)  given  by  (4.22), 


u> 

mx 

doe  g(oe)  =  1  -  f3  , 

foe 
mx 

1  3  1 

which  is  equal  to  7/8  for  f  =  -y  .  The  value  of  f  =  -j-  corresponds  to  con 
sidering  only  the  optical  modes.  (Half  of  the  modes  --  the  low-frequency, 
acoustical  ones  --  are  then  formally  neglected,  as  discussed  in  Sec.  V. ) 
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Note  that  feu  =  0. 79  ( 3. 8  x  10*3)  =  3. 02  x  10*3  rad/sec,  which  is  near  the 
13 

frequency  2 . 8  x  10  rad/sec  that  roughly  divides  the  optical  and  acoustical 
modes  in  Fig.  9. 

i  2  ^  _ 2 

Setting  f  =  y  in  (4. 25)-(4. 27)  and  using  a  =  a  -  a  and  u~/u:  =1.03 

^  u  j  mx 

for  NaCl  gives 

-  211  •  a,.  0.757  , 

a  2  =  0. 594,  aj  =  0.  594  -  (0. 757)2  =  0. 0210  =  (0. 145) 2  . 

3  1 

The  corresponding  values  for  f  =  -j-  are 

al  =  0.899  , 

a2  =  0.811  ,  OLy  =  0.0028  =  (0.053)2  . 

The  following  values  of  parameters  in  (4. 28)  are  used: 


m>  =  m^  =  5.  89  x  10-23g  , 
mr  =  2.32  x  10'23g  , 
nf  =  1.5, 

B  =  2.44  X  10**  dyne/cm2  , 


~  =  3. 82  X  10  g  , 

a  =2.  815  x  10  3 cm  , 
nn 

1/p  =  9.05  , 

e*  =  e  =  4. 80  x  10  *^*  esu  , 
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a),  =  3.09x  lO1^  rad/sec  ,  to  =  3. 8x  lO1^  rad/sec  ( /?  ’ s  independent 

1  of  u)  ) 

mx' 

oc  =  0.147  ,  Of,  =  0.757  , 

(Jl  1 

Dp  =  2.85  ,  K  =  5.92  x  103  cm"1  . 

7 

In  Fig.  1,  j39  •  •  •  0-  are  plotted  as  the  light  curves.  The  sum  0  =  L  0  ,  plotted 
11  n=2  n 

as  the  solid  heavy  curve,  shows  the  approximate  exponential  behavior  of  0 .  The 

agreement  between  this  theoretical  result  and  experimental  results  (points  and 
crosses  in  Fig.  1)  is  good  in  view  of  the  simplicity  of  the  model  and  the  assump¬ 
tions  made.  By  adjusting  the  values  of  two  parameters  in  (4. 28),  such  as  K  and 
D£ ,  the  experimental  data  can  be  fit  to  within  the  scatter  of  the  points.  By  changing 
the  value  of  the  single  interaction -strength  parameter  1/p  from  9.05  to  13.7,  the 
heavy  dashed  curve  in  Fig.  1  is  obtained.  Results  for  materials  other  than  NaCl, 
as  well  as  computer  evaluations  of  the  multiple  sums  of  0n  with  n  =  2  to  5  or  6, 
will  be  presented  in  a  future  report. 

Concerning  the  value  of  p,  it  should  be  mentioned  that  the  higher -order 
anharmonic  coefficients  are  much  more  sensitive  to  the  shape  of  the  potential 
curve  than  are  the  quadratic  terms,  which  are  used  in  the  analysis  of  the  binding 
energy  and  bulk  modulus.  Thus,  it  would  be  reasonable  to  replace  p  by  an 
n-dependent  factor  pn .  Since  there  appears  to  be  no  reasonable  method  of  de¬ 
termining  the  values  of  pn  from  first  principles  or  from  other  experiments  at 
present,  the  pn  would  have  to  be  considered  as  adjustable  parameters.  This 
would  give  five  parameters  to  adjust  in  order  to  fit  the  experimental  data  in 
Fig.  1.  Even  though  it  is  neither  surprising  nor  impressive  that  an  extremely 
good  fit  to  the  data  can  be  obtained  in  this  way,  it  is  interesting  that  the  value  of 
l/p„  is  equal  to  the  usual  value  of  1/p  and  the  values  of  1  /Py  1/p^  •  •  •  are  some¬ 
what  larger. 


if 
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The  n-phonon  regions,  marked  on  Fig.  1  for  n  =  3 - 6,  are  given  by  the 

relations  (n  -a )  <  to  <  ( n+ 1  -a )  ,  with  values  of  a  ranging  from  0. 5  for 

n  =  3  to  0.  8  for  n  =  6.  The  n-phonon  region  cannot  be  taken  as  ntOj  <  to  <  (n+  l)t0j , 

as  nte  <u)<(n  +  l)tc  ,  or  as  nto. -,  <  to  <  (n+ 1)  to.  „  ,  as  often  done  in  the 
mx  rax  LO  LO 

literature.  In  fact,  the  n-phonon  regions  shift  as  the  temperature  changes. 


It  has  not  been  established  conclusively  that  the  Lax-Bur  stein -Bom 
7  14a 

mechanism  is  negligible  in  polar  crystals.  The  fact  that  the  present 
anharmonic-potential  result  explains  the  observed  frequency  dependence  and 
order  of  magnitude  of  /3  suggests,  though  by  no  means  proves,  that  the  con¬ 
tribution  of  the  anharmonic-potential  mechanism  is  greater  than  that  of  the 
Lax-Burstein-Bom  mechanism  over  the  experimental  range  of  values  of  oo. 

The  absence  of  a  resonance  denominator  (2. 1)  in  the  Lax-Burstein-Bom 

-4 

result  would  eliminate  the  factor  to  in  (4. 28),  thus  giving  a  significant  devi¬ 
ation  from  an  exponential  frequency  dependence.  The  frequency  dependence 
and  numerical  values  of  the  matrix  elements  in  the  Lax-Bur  stein -Born  theory 
may  be  difficult  to  determine  accurately,  since  they  are  expected  to  be  rather 
strongly  model  dependent.  (A  simple  rigid-ion  model  gives  zero  for  the 
matrix  elements. ) 

1  3 

In  the  two  previous  reports,  *  numerical  values  of  /3  were  obtained 
by  making  approximations  that  were  different  from  those  of  the  present  report. 
Although  the  results  are  not  too  different,  the  present  results  are  believed  to 
be  the  more  accurate  of  the  two,  and  the  form  of  the  present  results  certainly 
is  more  convenient. 
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V.  The  Central-  Limit  Approximation 


The  central-limit  theorem  was  found  to  be  a  powerful  tool  for  obtaining 
the  values  of  0n  for  large  n ,  at  which  straightforward  computer  calculations 
are  impractical  at  present.  However,  considerable  care  must  be  exercised  in 
the  application  of  the  central-limit  approximation  lest  drastic  errors  should  be 
introduced.  In  the  present  section  and  in  Sec.  VI  it  will  be  shown:  that  the  con¬ 
fluence  processes  are  negligible;  that  including  the  confluence  processes  and 
using  the  central-limit  approximation  gives  drastically  incorrect  results  for  0 ; 
and  that  truncating  g(co)  can  give  improved  values  of  0  when  the  central-limit 
approximation  is  used  and  0(cc)  is  approximated  by  the  locus  of  the  peaks  of  the 
0n(cu)  curves. 


The  central-limit  approximation  will  be  applied  to  two  model  problems  in 
order  to  illustrate  these  difficulties  in  applying  the  approximation  to  the  present 
problem.  First,  consider  the  function  in  Fig.  10a .  The  nth  convolution  of 


f(oj)  =  fj  (a;)  is 


Vw> 


dtOj  f(co1) 


n 

f(w  )  6(w  -  E  a)  )  . 
n  i=l  n 


(5. 1) 


From  (5. 1)  it  is  easy  to  show  that 


y">  = 


fn<">  = 


dojj  f(co1)  f  (w  -Wj)  , 


dWl  fn-l  (aV  f  • 


(5.2) 


(5.  3) 
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The  result  of  the  successive  convolutions  in  (5. 2)  and  (5. 3)  is  shown  in  Fig.  10 
for  n  =  1,  2,  •  ■  •  5. 


The  central-limit  approximation  to  f  r(  co)  in  (5. 1)  is 


1 


f  Ju)  *- 

n  ^  2  IT  n  Mg 


exp 


(co  -  nco) 

2n  v 


(5.4) 


where 


co  =  l  d co  f(co)co  =  a  co  , 

mx 


(5.5) 


d co  f (co)  co2  =  a2co  2  , 
mx 


(5.  6) 


and 


2  2 
a!  =  OJ 


co2 


2  2 
a  co 
a  mx 


(5.7) 


Evaluating  the  integrals  (5.5)  and  (5.6)  for  f(co)  in  Fig.  10a  gives 


3 

TUU 

mx 


7  2 

<o  =  -4  oo 

12  mx 


“a  *  =  <°-“5"nuc>  ' 


The  result  for  n  =  5  is  shown  as  the  dashed  curve  in  Fig.  lOe .  The  agreement 

with  the  exact  solid  curve  is  poor,  since  the  central-limit  approximation  is  not 

valid  for  the  function  in  Fig.  10a  unless  n  is  very  large.  After  convolving 

f(co)  n  times,  the  Gaussian  half-width  of  each  peak  is  approximately  /rP  coo+, 

where  to  ,  =  co  N  48'  is  the  Gaussian  half -width  of  the  positive  half  of  f  (co) 
0+  mx 
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In  order  for  the  central -limit  approximation  to  be  valid,  the  width  J~\ T  W£y+ 
must  equal  the  spacing  (3/4)  of  the  two  sections  of  f(co),  roughly.  This 
gives  n/48  -  9/16,  or  n  as  27.  Thus,  n  must  be  quite  large  before  the  central- 
limit  result  is  valid  for  this  function  with  a  gap .  This  simple  example  will  be 
used  in  Sec.  VI  to  show  that  confluence  processes  are  negligible  in  the  expo¬ 
nential  frequency  region. 


Next  consider  the  function  F(a>)  in  Fig.  11a.  The  central-limit  approximation 
to  the  function 


F  (co) 
nv  ' 


F(u>j ) 


do) 

n 


F(U)  ) 
'  n ' 


n 

x  6  (  co  -  £  co.) 
i=l  1 


is 


F  (co) 
n'  ' 


1 


/ 2ffn  co„  exp 


( co  -  no?)2 

2na)(y2 


where  co  and  co  are  given  by  (5. 5)-(5.  7)  with  f(to)  replaced  by  F(o)).  The 
functions  p2(od)  and  F^^)  are  plotted  in  Fig.  11.  For  this  smooth  function 
with  no  "gaps"  the  central-limit  approximation  is  fairly  good  even  for  the 
very  small  values  of  n  =  2  and  3 . 

The  following  discussion  is  included  for  those  interested  in  the  details 

of  the  calculations.  As  a  rough  approximation  to  obtain  preliminary  results, 

the  curve  of  j3  =  L  8  as  a  function  of  co  was  approximated  by  the  locus  of 

i, 3 

the  peaks  of  the  individual  0n .  (The  size  of  the  error  introduced  is  apparent 
in  Fig.  1. )  Under  this  approximation,  it  is  found  that  formally  setting  a 
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Fig.  11.  Rectangular  function  F(to)  considered  as  the  sum  of  two  equal  rectangular 
functions  F^ou)  and  Fr(w).  The  dashed  curves  in  (b)  and  (c)  are  the  central-limit 
approximations  to  and  Fg ,  respectively. 
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low-frequency  portion  of  the  density  of  states  g  (to)  equal  to  zero  can  increase 
the  values  of  j3.  Increasing  the  value  of  f  in  (4.16),  (4. 19),  and  subsequent 
equations  corresponds  to  such  a  truncation  of  g(to).  The  reason  for  this  result 
that  decreasing  the  number  of  physical  processes  causes  an  increase  in  the 
value  of  |3  can  be  understood  as  follows.  In  Fig.  1  it  is  seen  that  the  value  of 
a  given  @n,  say  (3^ ,  at  its  peak  can  be  much  smaller  than  the  value  of  ^n_j  =  ft. 
at  the  same  frequency.  Truncating  g(a> )  causes  a  reduction  in  the  value  of  ft. 
(especially  at  its  peak  and  on  the  low-frequency  side  of  the  peak),  but  also  causes 
the  peak  to  shift  to  higher  frequencies.  This  new  peak  lies  nearer  the  true  curve 
of  0  =  Lpn  than  does  the  old  peak,  as  illustrated  by  the  following  example. 

In  order  to  examine  the  value  of  Fn(u,')  in  Fig.  11  in  the  wings  (away  from 
the  peak),  first  consider  F  (uj  )  as  the  sum  of  the  two  equal  rectangular  functions 

F(u>)  =  Fa(u>)  +  Fr(to)  ,  (5.8) 

as  shown  in  Fig.  11a.  The  first  convolution  F2(co)  is  then  the  sum  of  three  terms, 
as  seen  by  substituting  (5.8)  into  (5.2)  with  f(w)  replaced  by  F(a>).  This  gives 

Jco 

dWj  [F£(W1)  F^  (03  -  COj)  +  Fje(co1)  Fr(to  -Wj) 

+  Fr(co)  FA(w-co1)+  F^Wj)  Fr(w-w1)]  .  (5.9) 
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The  second  and  third  integrals  are  equal  as  seen  by  making  the  transformation 

f 

u  -  =  idj  in  either  integral.  Thus,  (5.  9)  gives 


F2(W)  =  F£2(CC)  +  2Fj *  Fr(a3)  +  Fr2(W)  * 


(5. 10) 


where  *  denotes  the  convolution.  In  the  convolution  notation,  F^^H  F^(W)  *  F^(UJ). 
The  three  functions  in  (5. 10)  are  plotted  in  Fig.  lib.  Similarly,  Fgfcd)  is  the  sum 
of  four  functions,  as  shown  in  Fig.  11c.  In  this  figure,  note  that  at  =  2.75, 

the  central-limit  approximation  (dashed  line)  has  a  smaller  value  than  the  exact 
result  (solid  line).  Also  note  that  at  this  value  of  00 /co^ ,  the  central-limit  ap¬ 
proximation  to  Fr2(t°)  is  a  closer  approximation  to  the  correct  value  than  is  the 
central-limit  approximation  to  F(co).  That  is,  if  F(cu)  is  truncated  to  Fr(co). 
the  resulting  central-limit  approximation  is  more  accurate  in  portions  of  the  wings 
than  the  central-limit  approximation  to  F  (co)  itself. 


In  the  same  way,  F  (<u)  could  be  truncated  to  the  rectangular  function  Frr(tc) 
extending  from  lc/(^mx  =  7/8  to  1 .  If  the  value  of  0n  in  the  wings  of  the  expo¬ 
nential  is  important,  formally  truncating  g  (to) ,  by  making  f  large  in  (4. 16)-(4. 19) 
for  example,  may  result  in  a  more  accurate  value  of  0n  .  This  result  was  the 
motivation  for  considering  two  values  of  f  in  Sec.  IV.  The  result  is  similar  to 
the  result  considered  above  and  in  Sec.  VI  that  the  central-liinit  approximations 
to  0  with  confluence  and  splitting  processes  included  is  not  as  good  an  approxima¬ 
tion  as  that  obtained  by  including  splitting  processes  only. 
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Sec.  B 


VI.  Effect  of  Confluence  Processes  on  B 

*  n 


It  will  now  be  shown  that  confluence  processes  are  negligible  in  the 
exponential  region.  The  analysis  of  Sec.  IV,  which  is  for  the  case  of  the 
splitting  of  the  fundamental  phonon  into  n  phonons  as  in  Fig.  3a,  can  be  easily 
extended  to  include  all  confluence  processes,  such  as  those  in  (b)  and  (c)  of 
Fig.  3.  The  general  confluence  process  has  m  output  phonons  and  n  -  m  + 1 
input  phonons.  The  method  below  sums  the  diagrams  over  all  values  of  m . 


It  can  be  shown  that  this  generalization  to  include  all  confluence  and  splitting 
processes  is  accomplished  simply  by  replacing  6  ( fl  -  u!q  )  by  6  ( ft  -0Oq  )  -  6  (SI  +  gCq ) 
in  (4. 6),  giving 


„  /  n  \  _  v 

ac±<°>  - 


2  n(ffl+l  lQ) 

r\  Wr\  - 


x  [  6  (  n  -  cOq  )  -  6  (  fi  +  coQ )  ] 


(6.1) 


and  defining  g(-co)  =  g(oo)  and 


n(-w)  -[n(w)  +  1]  ; 


(6.2) 


i.e.,  n(-uj)+l=  -n(co).  Standard  Green's  function  analyses  give  the  6  (J1  +  c0q) 
terms  automatically.  See  Eq.  (4. 19)  of  Ref.  6,  for  example.  The  resulting  expres¬ 


sions  for  (a  , )  are 
c± 


MQ> 


<ac±>  = 


<o  l  dO  Jj- 

mx  \  N 


^  n  (  fl)+ 1 


to. 


x  [  6(fi  -  <on)  -  6  (0+  con)  ] 


CO 


2n(cor.)+l 

mx  v,  Q  ^  /^v 

TT  £  - Z3T -  0  +  (Q) 

Q  Q 


(6.3) 


1 
3 


and  £  still  are  equal  to  the  n-fold  convolution  of  cr.  (O)  as  in  (4.5): 
ni  ± 


nc± 


~1T  [  d  «l  °c±  (  ftl> . j  dVc±  (fi„ 

“mx  4co 


)  6  (to  -  I  0*) 
i=l 


(6.4) 


Notice  that  the  summation  weighting  functions  a+(Ci)  in  Eq.  (4. 6)  of  Sec.  IV 
are  zero  for  negative  fi,  while  the  confluence  weighting  functions  in  (6. 1)  are  non¬ 
zero  for  both  negative  and  positive  values.  For  example,  if  the  angle  factors  0+(Q ) 
are  approximated  by  a  constant  $+ ,  then  (4. 6)  and  (4. 17)  give 

ot(0).  ^5L  f  d»g«0)  CO)  , 

a±  Jo 


6  CO 
±  mx 


o+(0)=  —r~~Z 
±  <oc  +  ) 


g(fl) 


n(  0)+l 

n 


for  >  0  , 


=  0  , 


'  Jg 


for  0  <  0  . 


(6.5) 


'lA'i 


acc.  d 


Similarly,  (6.1)  and  (4. 17)  give 


ac±  (ft)  = 


6$,  co  r 

*)  =  \ 
^ac±^  Jo 


dco  g(co) 


n  (Cl)  + 1 


co 


x  [  6(  Cl-w)  -  6(  C1+  u>)  ]  , 


°c± <n>  =  tztt  g(0) 

c±  ' 


for  Cl  >  0  , 


(6.6) 


_Pf  „|0|>  fcrO<0  • 


6$  .co 
±  mx 


Note  that  the  top  expression  in  ( 6.6)  is  formally  correct  for  Cl  <  0  also,  as  seen 
by  writing  Cl=  -  j  Cl  |  and  using  (6.2). 

The  central-limit  theorem  could  be  used  to  evaluate  the  multiple  integral 
in  (6. 4).  The  results,  however,  are  quite  misleading,  as  will  now  be  shown. 

It  is  necessary  to  carry  out  the  successive  convolutions  to  obtain  the  values 
°f  $2  »  ^3»  ’ ' '  •  Consider  £nc+  in  (6. 4).  In  the  example  above,  in  which 
0+(Q )  is  approximated  by  a  constant  $+,  a£+  ( Cl)  is  given  by  (6. 6).  In  the 
high -temperature  limit,  n((l)+l  -  n(Cl)  1/Cl;  thus,  if  g(Cl)  ~  Cl  over 
its  nonzero  region,  ac+(Cl)  is  constant  over  its  nonzero  regions.  The  function 
f(co)  in  Fig.  10a  is  therefore  a  reasonable  approximation  to  (Cl)  in  (6.6). 

In  fact,  the  ratio  of  1  to  3  of  the  heights  of  the  two  peaks  in  Fig.  10a  was  chosen 
with  foresight  to  agree  with  the  room -temperature  ratio  [n  (Cl)+l]/n(Cl)^  3  of 
the  two  peaks  in  on  the  positive  and  negative  Cl  axes,  according  to  (6. 6). 
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Physically,  the  positive-axis  peak  in  f(co)  in  Fig.  10a  corresponds  to  the 
creation  of  a  phonon,  and  the  negative-axis  peak  corresponds  to  the  annihilation 
of  a  phonon.  It  is  not  difficult  to  show,  by  considering  the  frequency  regions  over 
which  the  various  confluence  processes  can  conserve  energy,  for  example,  that 
this  interpretation  of  the  two  peaks  in  Fig.  10a  is  consistent  with  the  following 
interpretation  of  the  peaks  in  f  in  Fig.  10  b-e .  The  peaks  in  these  functions 
fn(co)  correspond  to  the  various  splitting  and  confluence  processes.  For  example, 
the  last  peak  on  the  right  of  f^  (u))  corresponds  to  the  splitting  of  the  fundamental 
mode  into  three  phonons  and  the  next  peak  to  the  left  corresponds  to  the  confluence 
of  a  fundamental  phonon  with  a  thermally  excited  phonon  to  form  two  phonons,  as 
illustrated  on  Fig.  10c . 

The  splitting  peak  (last  to  the  right)  is  not  the  largest  peak  for  a  given  n 

in  general.  However,  in  the  present  example,  the  other  peaks  can  be  neglected 

in  the  region  of  current  interest  since  they  are  obscured  by  the  contributions  of 

lower-order  terms  (smaller  n).  For  example,  the  last  (splitting)  peak  of  n  =  3 

gives  a  much  greater  contribution  to  /3  than  does  the  second  peak  from  the  right 

-1 

in  fcj  (to):  From  (2. 19)  with  the  frequency  dependence  of  (n^  + 1 )  neglected 
n 

and  £  -  {a  }  f _(co)/oo  ,  B  is  given  by  the  expression 

II  11  HJUv  11 

0n  ~  t  (<«>D  )"f  (co)  , 
n  to  n!  e  n 

where  (a)  —  0. 1 .  For  the  last  peak  in  Fig.  10c, 

<03>,  ~  Tf  l°'3  (0.633/^)  , 
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and  for  the  next-to-the-last  j3^  peak 

«>5 >2  ~  ^lo'5  (0.4285 /UnK>. 

Thus 

(ft.)  /(ft.)  *  3.38  x  10'4  «  1 
5  2  3  i 

so  that  (ft. )  is  negligible.  This  simple  example  clearly  illustrates  that  the 
confluence  processes  are  negligible  in  the  exponential  frequency  region;  only 
the  splitting  process  should  be  included. 

The  errors  introduced  by  formally  using  the  central-limit  approximation 
to  (6. 4),  which  includes  both  the  confluence  and  summation  processes,  are 
drastic.  Direct  calculation  shows  that  the  resulting  values  of  j3  are  incorrect 
by  orders  of  magnitude  in  some  cases.  The  difficulty  can  be  seen  simply  by 
considering  the  positions  of  the  peaks  of  the  various  0n  vs  u)  curves.  For 
example,  in  Fig.  lOe,  the  five-phonon  peak  is  at  co  as  2cof .  Experimentally, 
this  frequency  is  in  the  region  of  two-phonon  peaks,  not  five-phonon  peaks. 

An  even  more  dramatic  error  results  at  high  temperatures,  where  the  heights 
of  the  negative-axis  and  positive-axis  peaks  are  approximately  equal.  Then  all 
(for  reasonable  values  of  n  )  peak  near  u)  =  0 ,  which  is,  of  course,  below 
the  fundamental  resonance  frequency.  • 
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VII.  Effects  of  Terms  Other  Than  V^n+  ^  in  K 


In  the  previous  sections,  only  the  term  V^n+  ^  in  3C  [Eq.  (2.4)]  was 

considered.  In  determining  the  effect  of  the  other  terms  with  p  <  n  +  1 

in  (2. 4),  first  the  higher-order  terms  in  the  expansion  of  the  square  root  in 

R„  in  (2. 5)  are  needed.  With  the  definitions 
im  ' 


u  =  u.  /6.  , 

~  ~£m/  Am 

and,  formally, 

2 

u  •  u  =  u  , 

(2.5)  becomes 


Rlm  ^m  r ,  ,  2  ,  n  A  ,V2  , 

- - - =  [  1  +  u  +2  Au]  -1 

°m 


Expanding  the  square  root  gives 


^£m  uj£m  ’ 


6  •  u  =  Au  , 


R*m-6m  I1'2 


1 1/2  \  2  /  1/2  \  2  2  /  l/2\ 

=  ^  i  j(u  +2Au)  +  ^  2  )  <u  +2  Au)  +  |  j(u 


2  2  /  l/2  \  2  3 

i  +2Au)  +  ^  ^  j(u  +2Au)  +•• 

(7.1) 


where  the  binomial  coefficients  have  the  values 


1/2  \  , 

/  I/2  ] 

1  1 

jV2\ 

1 

1  I’  T 

\  2  1 

'  "  8  ’ 

(3  ) 

=  16  * 

l'2)  5 

I  1/2 

1  35 

I1'2] 

105 

4  /  128  ’ 

'  5 

1  1280  ’ 

(  6  ) 

5120 

(7.2) 


64 


Sec.  B 


Expanding  each  power  of  u  +2  Au  in  (7. 1)  and  collecting  equal  powers  of  u 
gives 

R.  -  6  ® 

im  v-  s 

— - ■  S  ru  .  (7.3) 

0  s=l 


where 


2A  =  A  , 


■  •  CD  C) 

-■no- -'•no' ■+ 
.■(?)(:)«*'•*  no 


(1-AZ)  , 


(2A)  =  -i- A (-1  + Az )  , 


=  -y  (-1  +  6 A2  -  5  A4)  , 


=  -i-  A  [3  -  10  A2  +  7  A4] 


r6  =  “IF  P  -  I5  A2  +  35  A4  -  21  A6  ] 


(7.4) 


Expanding  the  powers  of  (R.  -  6  ) / 6  gives  the  results  displayed  in  Table  I. 

jc  m 
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Multiplying  each  [(Rlm  -  6,m>/6imf  by  (-l)P(p!  pP)'1  and  summing 
gives 


K  =  C  L  iT1 

_ O  n 


(7.5) 


where 

C 


3 Ba  3p2 
nn  ^ 

1  -2  p 


and 


A  = 


2 !p2  2  I  p2 


1=- 


ll  +  2rlr2 


^3  = 


3  J  p 

4 

rl 


p22! 


4ip  3 !  p 


3!  p 

«  2  2 

3rl  r2  .  r2  +2rlr3 

X-  +  - - -  , 

“I 


(A3+  3pA3  -  3pA)  t 


2  1  p 


7" 1  4  [ a4+  3P  (5P+2)  A4  -  6p(l  +3p)  A2  +  3p2  1 
4 !  p  J 


V 


7  +  -V2  -  3(rl2f3+r22rl>  .  2<rlr4  +  r2r3> 

-5  '•  -4  77  3  - 5 - 

3  J  P  2  !  p 


5  !  p  4 !  p4 


—73  ['a5+  5p(2  -  9 p  -  2 1  p2  )  A5 
5 !  p 

-  10p(  1  -  6p  -  15p2  )  A3  -  15p2  ( 1  +  3p )  A  ]  . 
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%  -  TtV  frl6  -  30p  Tl\  +  »<»Vr3  +  6rlV+  3ii%> 

6:  p 

-  120p3(3r12r4  +  6r1r2r3  +  r23)  +  360p4  (2^  +  2  r2r4  +  r£  )] 


— l-z-  [ A6 ( 1+  1155p4  +  450p3  -  120p2  -  15p) 

6 !  p6 

+  A4(  -1875p4  -  675p3  +  240p2  +  15p )  A2  (765p4+270p3+ 180p2) 
+  (45p4  -  15p3)  .  r. 


The  first  term. 


,  n+1 


(n+l)I  p1^1 


in  each  in  (7. 6)  corresponds  to  the  value  used  in  the  previous  sections. 

It  is  tempting  to  approximate  by  1  //IP  in  (7.  6)  and  calcu¬ 

late  the  ratio  of  -^n  to  the  first  term  in  .  The  results  of  the  previous  sections 
could  then  be  multiplied  by  the  square  of  this  ratio  in  order  to  approximate  the 

effect  of  the  additional  terms.  Indeed,  this  increases  the  value  of  (|J  , ,)  /|3  ) 

n+1  mx  n  mx 

by  a  factor  of  approximately  3.7  for  small  n ,  giving  better  agreement  with  experi¬ 
ment.  However,  this  result  probably  is  not  correct  since  considerable  care  must  be 
exercised  in  approximations  of  this  type.  For  example,  if  A  is  approximated  by  1/3 

rather  than  by  l//1T,  the  value  of  (/3  . ,  )  /(jS  )  changes  erratically  with  in- 

mx  n  mx 

creasing  n. 
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From  *n  (7. 6)  it  is  seen  that  replacing  A3  by 
A3  +  3pA3  -  3p  A 

A 

takes  into  account  all  terms  in  3C„  fin  (2. 5)1.  Thus,  (  6_  •  u .  )' 

replaced  by 

(l  +  3p)(6  •  u .  )3-  3p(6  •  u.  )(u.  •  u.  ) 

v  K,v  m  ~jfcm  m  ~jtm  • 


The  corresponding  value  of  a 2  in  (2. 16)  is  replaced  by 


6 

=  L 
m=l 


(l+3p)(6m.  VVl'W 


V. 


-  »  <  «m  •  “f )  Bm«,>  •  um(Q2) 

For  uf  =  x,  (7. 8)  gives 

a2  =  Ux(Qi)Ux(Q2)-Pux(Qi)Ux(Q2) 

-pUx(Q1)Ux(Q2)  -  2 p  Uy(Q1)l^(Q2) 
-2pUz(Q1)Uz(Q2)  -  cc  . 


(7.7) 

in  (2.7)  is 


(7.8) 


(7.9) 
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Notice  that  the  coefficient  of  U  (Q.)U  (Q9)  from  the  extra  terms 

X  i  X  JL 
3 

[  3p  A  -  3p  A  in  (7.7)]  vanishes.  This  is  a  general  property  of  an  for  n  $ 
resulting  from  the  fact  that  all  =  0  when  A  =  1  for  n  2 .  See  (7. 4). 

For  a3,  (6m  •  U£m)4  n  (2.7)  is  replaced  by  [see  -$3  in  (7.6)] 

[l  +  3p(5p+2)](§m  •  U(m>4  -  6p(l +3p)(6m  •  ■  &lm) 

+  to2  '  W  • 

The  corresponding  value  of  a 3  in  (2. 16)  is 
6  (  ■ 

“3=  S,  ]  [l+3p(5p+2)](#m-  uf>  Um(Q  )  Um(Q2)  Um(Q  ) 
m=l  I 

-  3p(l+3p)  (6m  -  "uf )  Um(Q,)  [Bm(Q2)  •  Um(Q3)  ] 
-3p(l+3p,Um(Q1)Um(Q2)[S{  *  Hm<Q3>  1 

+3p2[uf  -  CHm<Q2)  -  Sim(Q3>  ]}  • 


(7.  10) 
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For  Uj  =  x ,  (7. 10)  gives 

a3  =  Ux(Ql>  Ux(Q3)  -3Pd+2p)[Ux(Q1)Ux(Q2)Ux(Q3) 

+  Ux(Q1)Ux(Q2)ax(Q3)  +  ^(QpUyCQ^UyCQg) 

+  uz<Qi>  uz(Q2}  +  3<°2  CUyCQi)  u£(Q2)  Uj(Q3) 

+  Uy(Qj)  Uy(Q2)  Uy  (Qg)+  Uz  (Qj)Uz  (Q2)  Uz  (Q  ) 

X  z  Z  X  X  X 

+  Uz(Q,)Uz(Q2)Uz(Q  )]  .  (7. 

x  *  y  **  y  ° 

These  results  (7. 9)  and  (7. 11)  are  being  used  in  the  computer  analyses, 
to  be  discussed  in  a  later  report. 
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VIII.  Higher-Order  Processes  In  The  n-Phonon  Region 


In  this  section  it  will  be  shown  for  NaCl  at  room  temperature  that  the 
higher-order  processes  shown  in  Figs.  2b,  2c,  12,  13,  14,  and  15  are  negligible 
for  reasonable  values  of  n.  The  results  suggest  that  higher-order  processes 
may  be  negligible  for  values  of  n  of  current  interest,  say  n  <  8 .  Additional 
analyses  will  be  made  elsewhere  to  establish  the  validity  of  this  suggestion. 


The  process  shown  in  Fig.  12a  will  be  called  the  single -vertex,  n-phonon 
process.  Processes  with  more  than  one  vertex  will  be  called  higher -order 
processes.  For  n  =  2  there  are  no  higher-order  processes  except  those  of 
the  type  in  Fig.  12b ,  which  can  be  accounted  for  simply  by  renormalizing  the 
phonons. 

For  n  =  3 ,  the  only  higher-order  process  of  interest  is  shown  in  Fig.  12c. 
The  contribution  to  frj  can  be  obtained  simply  as  follows:  By  standard  perturba¬ 
tion  theory,  the  matrix  element  for  the  process  in  Fig.  13a  (with  vertices 

aQ '  not  AQ  =  akb  +  a-kb  ^  1S 

'  %  /v  /v 


<fn  lKab+  li)=£4<1+2+3+f'  lA%)V 


t  t  I  ,+  V 

m23  am a2a3  1  1  m  f  > 


x  <lVY  |  A(Ka)  Vfmlaf.,  +  a+  |  l°2030m0f°>/fi«d  - -tum>  (8.1) 


where  |l-)E  |  n  ^  ±  1  >,  |nr)r  In  ),etc.,  and  permutations  of  2  and  3  and  of 
m  and  1  in  X  give  the  factor  of  4 .  Evaluating  the  sum  on  k  in  (8. 1)  by  using 
A(K&)  and  using  the  well  known  values  of  the  matrix  elements  gives 

0  (n  +1)V  00V.  . 

2  m  '  m2 3  fml 


<fn|KaW.!l> 


<fn|K3|i>  %  T  -«Jm)V. 


f  123 


(8.2) 
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n 


Fig.  12.  Phonon  relaxation  processes:  (a)  single-vertex,  n -phonon  process; 
(b)  processes  with  renormalized  phonon  lines;  (c)  three-phonon- summation 
higher -order  process. 
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where 


<  fn  |3C3  |  i  >  =  <l+2+3+f  |  A(Kfa)  Vfl23  |  l02°30f°  > 


The  diagram  in  Fig.  13b  gives  tne  same  result  as  (8.2),  but  with 
replaced  by  -com  and  nm  +  1  replaced  by  nm 


The  contribution  from  the  diagram  in  Fig.  13c  is 

< ^3+  I  1  *  =  E  2(nm+1)Vm23Vfml 
(fn  |3C3  1  i)  b  -SMw-tOj+u^JV^ 


(8.3) 


and  the  contribution  from  the  diagram  in  Fig.  13d  is  the  same  as  (8.  3),  but  with 

co  and  n  +1  replaced  by  -to  and  n  ,  respectively, 
m  m  m  m  J 


Adding  these  four  terms  gives 


<fa)J(3-22>1>  r  tomVfmlVm23 
<fn |Kj  |  i  >  bm 


The  value  of  this  ratio  in  the  three-phonon  region ,  say  at  to  =  3. 5  to  [  see  (5. 5)  ] 

is  of  interest.  With  this  value  of  to  and  with  to,  =  to  and  to  ^  to  ,  (8. 4)  re- 

1  m 

duces  to 


1^3-22  1 1 1  ^ _ 4 _  1 

<  fh  1 3C  |  i  >  15*75  hu2  b 


wmVfml  Vm23 


(8.5) 


The  value  of  =  vfQ  q  q  is  given  in  (2. 15).  This  result  (2. 15)  is 

12  3 

also  valid  for  V  ^  2  •  •  •  (n  + 1)  w*len  written  in  the  form: 
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Vl,2  •  •  •  n+1 


D  2 

3  p>/T  ,  nn  _ 1 _ 

~/T(l-2p)  ““  n!N*(n_1)  (wQ  ...OJ0  j171 

^1  Si+1 


x 


(8.6) 


where 


t 

a 

n 


6 

L 

p=l 


•  •  •  x 


U(rH-l)p 


(8.7) 


Substituting  (2. 15)  and  (8. 6)  into  (8. 5)  gives 


(  fn  Ik3_22  1 1  ^ 

<fn  |K3|"i) 


_6 _ 1_  Bann  /  1  \1/2  E  a2fmla2m23 

!°.5  l-2p  -2  ^m<m>  J  “3^23 


(8.8) 


1 


In  evaluating  the  sum  over  bm  in  (8.  8),  it  will  be  assumed  that  the  quasi  - 
selection  rule  that  the  matrix  element  for  the  splitting  of  a  fundamental  phonon 
into  two  optical  or  two  acoustical  phonons  in  NaCI- structure  crystals  is  negli¬ 
gibly  small.  Wavevector  conservation  of  the  three  output  phonons  will  be 
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neglected  formally,  and  the  value  of  (  fn  |  I  i)/(  fn  I  ^3  I  i)  will  be 


computed  at  co  =  3. 5  to,  assuming  that 


Wj  =  tc2  =  w3  =  OJ 


(8.9) 


for  the  three  output  phonons.  For  NaCl,  the  three  output  phonons  are  longi¬ 
tudinal  optical  (LO)  phonons.  Thus,  the  quasi -selection  rule  requires  that  the 
intermediate  phonon  m  be  an  acoustical  phonon.  The  case  of  m  being  a  longi¬ 
tudinal  acoustical  (LA)  phonon  will  be  considered.  Since  bm  =  LA  and  bj=  LO, 
a2fml  fr°m  (3*2)  can  be  approximated  by  (3.  9).  From  this  same  equation  (3.  9), 

with  6  •  u,  replaced  by  (m  / )'  U  to  account  for  the  fact  that  a„  00 
p  f  v  3  '  r'  <’  mp  2m23 

does  not  involve  the  fundamental  mode, 


/  mr  )>/2  I  _  /“< 

a2m23  “  \  m<  /  l  u<krnLAx  " 


m^  U>k  LAx  6 
^  m 


ik_„  a 
mx  nn 


X  ^^l^LOx  J 

X  (u<k3LOx  "  J  m~.  u 


ik2xann 


>k2LOx 


>kg  LOx 


ik„  a  \ 

3x  nn  \ 

e  1  -  cc  . 


(8. 10) 


As  in  Sec.  Ill,  we  set  u<LOx  =  1/''r3\  u>LOx  =  °*  u<LAx  =  °’  and 
u>LAx  =  •  These  values,  with  (8. 10)  and  sin2^^  =  y  formally. 


2  mr 


2m23 


(8.11) 
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Next  consider  the  value  of  a  3^23 
indicate  that 


The  same  arguments  that  gave  (8. 11) 


(8. 12) 


is  a  reasonable  approximation.  Substituting  (3.  9),  (8. 11),  and  (8. 12)  into  (8.  8) 
gives 


(  fa  1  K3_22  l ^ 

LOLA 

/  „  .2  2 

I  /  1  Bann  \  m<  mr 

<fn 

*3 

i> 

[20.7^=^^  m>3  j 

=  1. 0  x  10'2  .  (8. 13) 

Even  if  the  other  two  acoustical  intermediate  states  contribute  as  much  as  the 

LA  state,  the  total  ratio  in  (8. 13)  is  increased  only  by  a  factor  of  3,  giving 
-2 

3  x  10  for  the  value  of  this  ratio.  Thus,  the  order -of- magnitude  calcula¬ 
tion  suggests  that  the  higher -order  processes  are  negligible  in  NaCl  for  n  =  3. 

In  considering  larger  values  of  n ,  several  classes  of  diagrams  can  be 

studied.  First  consider  the  n-phonon  diagram  in  Fig.  14a.  The  energy  de- 

2  2 

nominator  corresponding  to  (to  -u^  )  -  tom  in  (8.4)  is 
. %-2)2-“m2 

Assuming  that  w  =  (n  +  -j)  £0  and  Wj  =  co2  =  •  •  •  =  wn_  j  =  tom  =  to  ,  this  energy 
denominator  reduces  to 

1  —  _  2  _  9  _o 

[  n+ yW  -(n-2)  «]  -  to  =  5.25  to  . 
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The  rest  of  the  analysis  is  the  same  as  that  leading  to  (8. 8),  and  the  results 
are  the  same  as  (8.  8)  with  the  appropriate  changes  in  the  matrix  elements: 

(  fa  l^n-2(n-l)  [  *  ^  _  _ 6_  /_n_\  1  ^ann  /  1  \  —  tt(n-l)fml‘  •  ,a2mnn-l 

10. 5  \  3  /  1  - 2p  -2  l  m<m>  )  b  “nfl2*--n 

'  '  m 

(8. 14) 


The  value  of  .  is  given  by  (8. 11): 

L  mn  n-1  ° 

2  2  mr 

|a2mnn-l  I  _  27  m> 


(8. 15) 


With  all  output  phonons  1  •  ■  •  n  as  LO  phonons  and  the  intermediate  phonon  as  a 
LA  phonon,  the  arguments  leading  to  (3. 9)  indicate  that 


a 


nfl2 


■n 


2 

I  a(n-l)fml  •  *  •  | 


(8. 16) 


are  reasonable  approximations.  A  factor  of  2/3  was  included  in  (8. 16)  as  dis¬ 
cussed  in  Sec.  III.  From  (8. 14)-(8. 16)  with  £  =  3  formally. 


<fn  |K 


~n-2(n-l) 


|i> 


<  fn  |3Cn  |  i) 


(8. 17) 


2 

This  result  (8. 17)  is  (n/3)  times  larger  than  the  n  =  3  resu’t  (8. 13).  The  ratio 
in  (8. 17)  is  equal  to  unity  for  n  =  190.  Thus,  it  is  clear  that  the  contribution  to  ft 
from  the  diagram  in  Fig.  14a  is  negligible. 
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Next  consider  the  diagram  in  Fig.  14b.  For  n  =  3,  this  diagram  reduces  to 
that  of  Fig.  12c,  and  the  matrix  element  is  given  by  (8. 13).  For  larger  values 
of  n ,  a  straightforward  evaluation  of  the  sums  over  the  diagrams  corresponding 
to  those  in  Fig.  13  would  be  tedious.  A  simple  method  of  evaluating  the  sums  is 
first  illustrated  for  n  =  3 .  Since  the  ratio  (  fn  1 _22|i>/<fe|3<3|i>  is  inde¬ 
pendent  of  temperature,  ail  occupation  numbers,  such  as  nm  in  (8. 2),  must  sum 
to  zero.  Thus,  only  terms  containing  factors  of  (n+ 1)  need  be  considered.  For 
example,  the  diagrams  in  Figs.  13(a)  and  (c)  are  considered,  and  those  in  Figs. 
13(b)  and  (d)  are  dropped.  [  For  the  m  arrow  going  from  left  to  right,  the  term 
aman^  gives  rise  to  nm+  1 ,  while  the  m  arrow  going  from  right  to  left  cor¬ 
responds  to  a  =  n  .1 
1  mm  m  J 

Since  the  sum  over  the  four  diagrams  contains  the  factor  co ,  the  energy 
denominators  can  be  expanded,  and  only  the  com  terms  retained.  For  example, 
the  energy  denominator  corresponding  to  Fig.  13(a)  is 


1  1  ,  m 

“  -  “l  -  %.  ■  “i  (^7 


(W  -  COj  ) 


(8. 18) 


and  that  for  Fig.  13(c)  is 


<“-“x+“n 


CO  -  CO. 


1  (W-Wjf 


(CO  -  COj  ) 


(8. 19) 
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For  large  n ,  n  +  -j-  a:  n .  Thus  we  set  to  =  (n  +  y-  )  to  ^  n  to ,  or  to  =  3  to  for 
n  =  3  in  anticipation  of  the  large  values  of  n  of  interest.  Adding  (8. 18)  and 
(8. 19)  and  using  to  —  3  to  and  tOj  —  to  gives 


2—2 
(3.5  -  If 


(8.20) 


which  corresponds  to 


(to  -  cOj )  -ton 


(8.21) 


in  (8.4).  For  large  n,  the  factor  tom  in  the  denominator  of  (8.21)  is  negligible, 
and  (8.20)  and  (8. 21)  are  equivalent.  To  shorten  the  notation  further,  set 
to  =  3  tOg,  tOj  =  00^  and  tom  =  mto^  in  (8. 18)  and  (8. 19).  This  reduces  (8. 18) 
and  (8. 19)  to 


1 

2  -m 


1  m  , 

^T+W  °  T+ 


and  the  sum  is  to  2m/4 ,  which  can  be  written  as 


2m  _  2m 

(2I)2 


(8.22) 


For  n  =  4 ,  a  straightforward  evaluation  of  the  sum  over  the  energy  denomi¬ 
nators  of  the  24  diagrams  gives 


to  to 

mj  m2  mj rru  (2)  mjm^ 


(8.23) 
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This  result  (8. 23)  can  be  obtained  by  the  method  above  as  follows:  The  energy 
denominators  corresponding  to  the  diagrams  in  Figs.  15(a)-(f)  are,  respectively 


Sec.  B 


Thus 


^  I  Kn-22  •  •  •  1 1  ^  2n  '  2 


=  + 


<  fit  1 3Cn  |  i  >  [(n-l)!]2  (hw2)n'2  bm  b, 


w  •  •  •  <0_ 


ml  “mn-2  1 


m,  m. 


n-2 


V-.  Vm  m2***Vm  /„  1\H 

flmj  nij m£  mn_2(n-l) 

- 

fl2  •  **n 


(8.25) 


Substituting  (2. 15)  and  (8.6)  into  (8.25)  gives 


<ft|*n-22-  i‘>  . 

<  fa  |  3C3  |  i ) 


_ n _ 

2(n-  1)! 


Ba 


nn 


2(1-2 p) 


n-2 


a2flm10£2m.mo2  ",0£2m  0(n-l)n 
x  E  •••  E  _ - _ — _ — _ 

bml  bm(n-2)  anfl2  *  *  *  n 


(8.26) 


For  n  =  3,  the  numerical  coefficient  9/8  in  (8.26)  is  slightly  different  from  the 
corresponding  value  9/10. 5  from  (8.8).  The  difference  is  that  tom2  in  (8.21) 
was  not  neglected  in  the  derivation  of  (2. 15)  and  n  +  y  was  not  approximated 
by  n. 


The  values  of  the  multiple  sum  in  (8. 26)  can  be  estimated  by  formally 
setting 


2 
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and  using  (8. 11)  and  (8. 16),  which  gives 


S  ^  lK3-22‘*  •  1  ^ 

~  n  / 

(  5/21  Bann  \ 

i  <  fn  |3C3  |  i  > 

1 

to 

4k 

13 

*-• 

■ — ' 

2(1  -2p)  1 

2n-4 


(8.27) 


Since  the  value  of  the  ratio  in  (8. 27)  decreases  with  increasing  n  and  is  negli¬ 
gible  for  n  =  3 ,  the  order  of  magnitude  calculation  suggests  that  the  contribution 
to  ft  from  the  diagrams  in  Fig.  14b  is  negligible. 

Although  the  temperature  and  frequency  dependence  of  j8  for  the  higher - 
order  processes  is  of  academic  interest  only  since  these  processes  are  negli¬ 
gible,  it  is  mentioned  that  the  frequency  dependence  of  the  higher -order  pro¬ 
cesses  is  different  from  that  of  the  single -vertex  process.  The  formal  expression 
for  the  temperature  dependence  of  the  higher -order  phonon  processes  is  the  same 
as  that  of  the  single-vertex  n-phonon  process.  However,  wave-vector  conserva¬ 
tion  in  the  intermediate  states  may  cause  the  final -state  phonons  to  be  different 
in  the  single-vertex  and  higher-order  processes,  thus  giving  some  difference  in 
the  temperature  dependence  of  the  two  types  of  processes.  Furthermore,  the 
factors  B,  a^ ,  and  co  in  (8.27)  are  temperature  dependent.  At  the  present  state 
of  knowledge  ( of  the  phonon  dispersion  relations  and  the  coefficients  in  the  scat¬ 
tering  Hamiltonian,  for  example )  these  temperature  differences  could  not  be 
used  to  distinguish  between  the  two  types  of  processes. 
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IX.  Anharmonic -Einstein -Oscillator  Model 

A  study  of  the  absorption  of  electromagnetic  radiation  by  a  collection  of 
non-interacting  but  anharmonic  oscillators  has  been  initiated  as  an  alternate 
approach  to  the  multiphonon  absorption  problem.  A  ^eliminary  report  of  the 
results  is  presented  here.  A  complete  presentation  will  be  given  in  a  future 
report. 

Since  most  of  the  experimental  measurements  of  interest  to  this  study 
have  been  carried  out  at  room  temperature  or  higher,  purely  classical  methods 
were  used.  The  first  step  was  to  derive  an  expression  for  the  absorption  coef¬ 
ficient  of  the  oscillator  array,  assuming  the  system  in  thermal  equilibrium  at 
temperature  T .  The  resulting  general  form  may  be  applied  easily  and  quickly 
to  compute  the  absorption  coefficients  for  a  number  of  anharmonic  potentials  of 
physical  interest. 

The  Liouville  equation  for  the  classical  distribution  function  p(q,p,t)  of 
the  oscillator  array  in  the  presence  of  a  spatially  uniform  driving  field  of  fre¬ 
quency  a)  was  treated  by  perturbation  theory  to  obtain  a  classical  form  of  the 
well  known  Kubo  formula  for  the  absorption  coefficient.  The  Kubo  formula  was 
then  reduced  to  a  form  that  may  readily  be  computed  once  the  nature  of  the  par¬ 
ticle  orbit  in  the  anharmonic  potential  is  known.  The  final  form  is  particularly 
convenient  for  use  in  the  quasi-harmonic  region,  where  anharmonic  corrections 
to  the  particle  orbit  are  small.  This  is  the  temperature  region  of  interest  in  the 
discussion  of  both  solids  below  their  melting  temperature,  and  for  molecular 
gases  below  the  temperature  at  v  hith  dissociation  occurs. 
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The  formula  was  applied  to  the  case  of  particles  whose  motion  is  governed 

by  the  Morse  potential,  an  anharmonic  potential  used  frequently  in  molecular 

physics.  In  this  quasi-harmonic  region,  we  obtain  a  particularly  simple  result 

for  the  absorption  coefficient.  If  the  fundamental  vibration  frequency  of  the 

th 

molecule  is  then  the  intensity  of  the  n  harmonic  absorption  at  OJn  =  n 
is  given  by  the  simple  expression 


where  D  is  the  dissociation  energy,  and  Pj  the  strength  of  the  fundamental  ab¬ 
sorption  line  at  u>q.  For  temperatures  near  room  temperature,  kT  «  D, 
and  a  plot  of  Pn  against  n  gives  a  form  very  close  to  the  exponentially  decaying 
form  observed  in  many  crystalline  materials. 

Mangir  and  Hellwarth^  have  treated  the  same  problem  from  a  quantum 
mechanical  point  of  view  and  obtained  a  result  that  reduces  to  ours  in  the 
correspondence-principle  limit  fi  -*  0 .  A  detailed  comparison  between  our 
calculation  and  theirs  should  prove  most  illuminating. 

We  believe  a  form  similar  to  the  result  quoted  above  will  very  likely 
emerge  for  large  n  from  the  study  of  any  potential  for  which  the  deviations  of 
the  small-amplitude  motions  from  the  harmonic  approximation  are  small.  It  is 
interesting  to  contrast  this  near -exponential  behavior  with  that  appropriate  to  a 
potential  for  which  the  harmonic  approximation  cannot  be  applied.  An  example 
of  such  a  potential  would  be  a  simple  square-well  potential,  with  a  flat  bottom 
and  infinitely  steep  sides.  Applying  the  method  to  this  problem  indicates  that 
at  high  frequencies  the  absorption  does  not  fall  off  in  a  near  exponential  manner, 
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2 

but  rather  as  1/co  .  We  plan  to  examine  two  or  three  more  potential  functions 
to  see  whether  or  not  one  can  make  a  clear  systematic  statement  about  the  rela¬ 
tion  of  the  form  of  the  potential  function,  and  the  frequency  dependence  of  the 
absorption  coefficient  at  high  frequencies. 

The  comparison  between  the  results  for  the  Morse  potential  and  the  square 
well  leads  us  to  conjecture  about  the  possible  role  of  certain  impurities  in  in¬ 
creasing  the  linear  absorption  coefficient  of  crystals  above  die  value  expected 
for  the  intrinsic  processes  apparently  responsible  for  the  near -exponential  tail 
at  high  frequencies.  It  is  well  known  that  in  the  alkali  halides,  certain  impuri¬ 
ties  (Li  and  Ag,  as  examples)  move  in  potential  wells  that  are  very  flat,  since 
very  low,  strongly  temperature -dependent  frequencies  are  associated  with  these 

impurities.  It  may  be  that  the  presence  of  these  impunities  will  give  contribu- 

_2 

tions  to  the  absorption  coefficient  which  falls  off  like  to  rather  than  exponen¬ 
tially  with  frequency.  Consequently,  at  high  frequencies,  the  value  of  the 
absorption  coefficient  may  be  greatly  affected  by  these  impurities. 

The  results  obtained  to  date  will  enable  us  to  make  semi-quantitative 
estimates  of  the  impurity  concentrations  required  to  give  rise  to  a  given  value 
of  the  absorption  coefficient.  The  simplicity  of  the  method  may  make  it  useful 
in  other  problems,  such  as  examining  nonlinear  absorption  mechanisms. 
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C.  THEORY  OF  INFRARED  ABSORPTION  AND  MATERIAL 
FAILURE  IN  CRYSTALS  CONTAINING  INCLUSIONS 

Small  amounts  of  macroscopic  inclusions  in  or  on  the  surface  of  a 

-7 

nonabsorbing  host  crystal,  in  some  cases  a  fractional  volume  of  ~  10  to 

10  can  give  rise  to  a  value  of  the  optical  absorption  coefficient  jS  of  10  * 

cm  1  (a  typical  value  of  current  interest).  For  various  types  of  inclusions,  the 

frequency  dependence  of  /3  ranges  from  increasing  as  to  »  to  independent  of  to, 

to  exponentially  decreasing  with  to.  The  temperature  dependence  ranges  from 

independent  of  T  to  increasing  as  T*5  in  the  high -temperature  limit,  where 

p  =  2  -  4  typically.  Simple  expressions  for  the  absorption  cross  section  are 

derived  for  various  cases  of  practical  interest.  The  cross  sections  are  used 

to  derive  expressions  for  0  for  the  four  cases  of  large  inclusions  of  strong 

and  weak  absorbers  and  of  small  inclusions  of  dielectric  and  metallic  particles. 

A  result  of  general  interest  in  the  theory  of  scattering  is  that  the  absorption 

2 

cross  section  of  large  spheres  is  equal  to  (1  -  R^  )rra  ,  where  Rj_  is  the  re- 

flection  coefficient  for  normal  incidence,  rather  than  (1  -  (R))na  ,  where 

(R  )  is  the  geometric -opucs  average  reflection  over  the  illuminated  hemisphere. 

The  material  failure  resulting  from  local  heating  of  macroscopic  inclusions  is  a 

far  greater  problem  in  high-intensity  short-pulse  systems  than  in  low-intensity 

long-pulse  or  cw  systems  having  equal  average  intensity.  Microsecond  pulses 

2 

with  energy  densities  as  low  as  ~5  J/cm  can  cause  material  failure. 
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I.  Introduction 

Several  theoretical  studies  of  intrinsic  infrared  absorption  are  currently 
underway. 1 5 '  *  *  1  >  ®  it  is  expected  that  the  results  of  these  studies,  especially 
in  conjunction  with  experimental  investigations,  will  afford  estimates  of  the 
values  of  /3  for  various  materials  that  could  be  approached  by  removing  the 
extrinsic  sources  of  absorption.  The  present  investigation  is  concerned  with 
another  aspect  of  the  problem  --  that  of  extrinsic  absorption  by  macroscopic 
inclusions.  Since  the  inclusions  may  be  either  in  the  bulk  of  the  crystals  or 
on  their  surfaces,  the  results  are  of  interest  in  surface  studies  as  well  as 
studies  of  bulk  properties. 

Absorption  is  of  greater  interest  than  scattering  in  the  study  of  high-power 

-4  -1 

infrared  windows.  A  value  of  ^  =  10  cm  for  the  absorption  coefficient  can 

cause  sufficient  heating  for  the  system  to  fail,  for  example,  by  thermal  defocus- 

18*20  *4 

ing  or  by  window  fracture.  But  an  equal  amount  of  scattering  0gcat  =10 

cm  1  may  be  tolerable.  Thus,  scattering  will  be  neglected  here,  except  to 

mention  that  the  observation  of  the  associated  scattering  may  help  to  identify 

21 

the  absorption  mechanism  in  some  cases.  Winsor  has  shown  theoretically 
that  scattering,  especially  in  conjunction  with  total  internal  reflection  at  the 
host  crystal  boundaries,  may  increase  the  measured  value  of  j3  by  increasing 
the  path  lengths  of  the  rays  in  the  crystal. 

22 

There  is  great  interest  in  lowering  the  value  of  jS  for  candidate  ma¬ 
terials  for  high-power  infrared  windows.  Materials  of  current  interest  include 
ZnSe,  CdTe,  KC1,  KBr,  and  TI 1173  glass  (Ge2gSbj2  Se^).  Results,  such  as 


/ 
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those  presented  here,  concerning  the  magnitude  and  the  temperature  and 
frequency  dependence  of  /3  for  extrinsic  and  intrinsic  processes  are  needed 
in  the  material-improvement  programs.  The  effects  of  inclusions  on  absorp¬ 
tion  in  the  infrared  have  not  been  studied  experimentally.  Such  experimental 
studies  would  be  welcomed. 

It  has  been  realized  for  some  time  that  windows  in  high-power  infrared 

laser  systems  can  cause  the  system  to  fail  with  or  without  actual  mechanical 

18-20 

failure  of  the  window.  Examples  are  thermally  induced  fracture  of  the 
window  material  and  thermal  defocusing  of  the  laser  beam  by  the  window 
resulting  from  uneven  heating  of  the  window,  respectively.  The  increase  in 
the  spatial  average  of  /3  resulting  from  inclusions  can  cause  either  type  of 
failure.  In  addition  to  this  overall  heating  of  the  material,  localized  heating 
in  the  vicinity  of  inclusions  can  cause  material  failure.  This  local -heating 
failure  effect  is  considered  in  Sec.  V. 

For  sufficiently  large  concentrations  of  inclusions,  the  localized  heating 
can  also  cause  considerable  optical  degradation  of  the  beam.  Li  the  present 
study  it  is  assumed  that  the  concentration  of  inclusions  is  so  small  that  this 
localized-heating  type  of  optical  distortion  is  negligible.  However,  it  should 
be  mentioned  that  scattering  may  be  considerably  larger  at  high -power  levels 
than  at  low  levels  as  a  result  of  the  scattering  by  the  heated  host  material  near 
the  inclusion.  This  effect  should  be  greatest  for  scattering  near  the  forward 
angle. 
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D.  Analysis  of  Absorption  Cross  Sections 


Consider  a  sample  consisting  of  a  nonabsorbing  host  material  of  dielectric 
constant  (real)  containing  Nj  inclusions  per  unit  volume,  each  inclusion 
having  absorption  cross  section  0  .  •  Multiple  scattering  will  be  neglected  -- 
a  reasonable  approximation  for  the  present  case  of  small  impurity  concentrations. 
The  well  known  result  for  the  absorption  coefficient  0  of  the  sample  is  then 


*  =  "abs  NI 


(2.1) 


Only  the  case  of  spherical  inclusions  of  radius  a  will  be  considered.  The 

absorption  cross  section  is  not  generally  equal  to  the  geometrical  cross  section 
2 

it  a  .  For  ka  «  1 ,  where  k  =  2  ff/A  ^  ,  with  A^  the  wavelength  of  the  radi- 

2 

ation  in  the  host  material,  the  value  of  typically  is  small  (CTabg  «  it  a  ). 

In  the  case  of  ka  »  1  and  |  €  |  »  1 ,  where  c  =  €j/ch,  with  €j  the  dielectric 

constant  of  the  inclusion,  the  reflection  at  the  surface  of  the  sphere  is  great, 

2 

which  again  makes  «  IT  a  . 


Simple  expressions  for  for  various  limiting  cases  will  be  derived 

23,24 

from  the  classic  result  of  Mie: 

°abs =  -jr  (2£  +1)  [*e  (a*  +  V  “  |  al  |  "  |  bl  |  J  '  (2' 


2) 


where  Re  denotes  "the  real  part  of . "  The  values  of  a^  and  b^  in  the  regions 
of  interest  are  given  below  in  Eqs.  (2. 3),  (2. 7),  and  (2. 8). 
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24,25 

For  small  spheres  (ka  «  1 ),  the  sums  are  well  approximated  by 
the  single  term  containing  a^  ,  where 


al  "  -i  3  €+2  * 


The  time  dependence  is  taken  as  exp(-ioot).  Some  authors  use  exp(+ia>t),  in 
v'hich  case  all  factors  of  i  in  (2. 3)  and  below  would  be  changed  to  -i .  With 

i  |2 

e  =  +  i  ,  substituting  (2. 3)  and  |  aj  j  «  Re  a j  into  (2.2)  gives 


12  ka 


<v2)  +  v 


tra2  , 


for  ka  «  1 


Two  limiting  cases  of  (2.4)  are  of  interest.  For  ,  which  is 

satisfied  for  nonmetals  at  frequencies  not  too  near  the  fundamental  resonance 
frequency  or  the  high-frequency  absorption  edge,  and  Cj  are  related 
the  expression 


-  2 c9m(n)  k  ^  R__ 


(2.5) 


Using  this  result  to  eliminate  Cjk  in  (2. 4)  gives 


.  12  Re  n  /a  v  _  2 

*abs  “  ’  T  * 

s  (ch+2  r 


€j  <<:  €ft  ‘ 


The  second  limiting  case  of  (2. 4)  is  that  of  small  metallic  inclusions. 
The  Drude  expression  for  the  dielectric  constant  is 


€  =  €_  - 


€,,(60  +  iCoT  ) 


(2.7) 
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where  T  is  the  electron  relaxation  frequency  (often  written  as  1/T),  is 

2  1/2 

the  contribution  to  €  from  the  core  electrons,  and  =  (4irNe  /m)  '  is  the 

plasma  frequency,  which  has  a  typical  value  of  =  5  x  lO1^  rad/sec.  At 

14 

10. 6 Jim,  oj  —  1.9  x  10  rad/sec.  There  are  two  contributions  to  the  relaxa¬ 
tion  frequency  T:  T  =  +  rSu-  where  the  bulk  contribution  has  a  typical 

value1  ^  of  I'BUa5x  10  ^  sec  1.  The  value  of  the  surface -scattering  contribu- 

Q 

tion  Tgu  is  ~Vp/a ,  where  the  Fermi  velocity  Vp  has  the  value  vp  a  10 

13  - 1  g 

cm/sec  for  many  metals.  With  rB„*5*  10  sec  and  vp  s*  10  cm/sec, 

Tc  >  for  a  ^  200  A . 

Su  Bu 

2 

For  typical  metals  at  10. 6 Jim,  »  |  co  +  icoT  | ,  and  (2.7)  gives 
2  2 

“p  wp  /  ,r\ 

•  ^  CH  —  ~  2  ~  ~  2  "  5”  (  ^  ^  qj  )  »  ^  c  ^  ^  ~  ^  a  • 

H  w^+itor  ur  +  r*  '  w  R  R 


Both  the  real  and  imaginary  parts  of  the  dielectric  constant  are  large  in  magni¬ 
tude  for  small  particles  of  typical  metals  at  10. 6  Jim .  Substituting  (2. 8)  into 
(2.4)  gives 

3/2  2  p  o  9  9 

ffabs  =  12  fH  ^7  ™  ■  “p  ^  1“  +i“r  I  •  &i> 


for  metals.  This  simple  expression  for  the  absorption  cross  section  will  be 
used  in  Sec.  £Q. 


Next  consider  the  case  of  lars 


3,  so  that 


ka  »  1  , 


T  PI 


0Ta  »  1 


(2. 10) 
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are  satisfied,  where  =  2  k  Jm  n  is  the  bulk  absorption  coefficient  of  the 
inclusion,  n  =  /T1 ,  and  denotes  "the  imaginary  part  of."  Ihe  inequality 
j  »  1  typically  is  satisfied  for  metals  and  strongly  absorbing  dielectrics. 
It  will  be  shown  that 


abs 


=  (1  -  R±)ira  , 


ka,  y  fya  »  1 


(2.11) 


where  =  |  (n  -  1 ) / (n+ 1 )  |  is  the  reflection  coefficient  for  normally  inci¬ 
dent  radiation. 


The  asymptotic  expressions  for  the  a^  and  b^  in  (2.2)  are 


25a 


a„  - 


.  I  -i  ka 
-l  e 


sinx  -ncosx  tany 
1  -  in  tany 


25  0, 


for  i  ^  ka 
for  l  >  ka 


b,  =  it+1e-ika  fori  Ska 

*>  tany  +  in 


*  0, 


for  i  >  ka 


(2. 12) 


(2. 13) 


where  x  -  ka  -  y  (i  + 1 )  jr  and  y=nka  -  y  (£  + 1 )  fl  .  The  approximation 
tany  —  i  is  valid  when  «9m(n)ka  »  1  is  satisfied,  as  seen  from  the  expres¬ 


sion 


Kn(A+1B>  -  WrfrSHA  *  lta"h  2B  a  ‘  ’ 


(2. 14) 


where  the  approximations  are  valid  for  B  »  1 . 


Using  tany  -  i  and 


,j£+l-ika 
i  e 


JL+1 

-  (•**) 


-ika  -ix 
e  -  e 


(2. 15) 
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reduces  (2.12)  and  (2.13)  to 

_  ..  _-ix  i  sinx  +  ncosx 

h  - e  — 1+5 — 

for  1  ^  ka .  From  (2. 16),  it  is  easy  to  show  that 

aJ2  +  I  |2  *  -  al  +  bi2:1-  <2-17) 


br 


,-ix  cosx+insinx 
1+n 


(2. 16) 


for  1  *  ka,  and  from  (2.12)  and  (2.13), a^  a  2:  0  for  1  >  ka. 


Substituting  (2. 17)  into  (2. 2)  with  the  upper  limit  in  the  sum  replaced  by  ka  and 
using 


1 


i  i2 

I  n  1+1  _  2  Re  n  _  1 
|n+l|2  I  n+ 1 1 2  2 


(1  -Rj^) 


gives 


abs 


^  2tr  1 

”  7  T 


ka 

(1-R.  )  L 

-L  1=1 


(21+1)  . 


(2. 18) 


(2.19) 


This  expression,  with 


ka  2 

£  (2X  +  1 )  a:  (ka)  ,  for  ka  »  1  ,  (2.20) 

1=1 


gives  (2. 12). 

When  -j-  Hja  >>  1  is  not  satisfied,  that  is,  when  a  ray  traversing  the  inclu¬ 
sion  is  not  completely  absorbed,  it  is  reasonable  that  (2.11)  should  be  replaced  by 
the  approximate  expression 

aabs  S  U-R1)(l-e'^ia)ira2  .  (2.21) 
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The  well  known  result  that  the  extinction  cross  section  <x  =  +  oscat 

is  equal  to  twice  the  geometrical  cross  section  when  ka  »  1  follows  directly 
from  (2. 17).  Substituting  =  1  for  l  ^  ka  and  =  0  for  l  ^  ka 

into  the  expression 

"ex,  *  ^  jj  <“+I>Re(“i  +  bi) 

2 

and  using  (2..20)  gives  or^  =  2  ft  a  . 

For  large  spheres  with  (2. 10)  satisfied,  it  is  not  difficult  to  show  that  using 
geometrical  optics  approximations  gives  the  result 

aabs  =  (1  '  <R>)ffa2  »  (2.22) 


where  (  R  )  is  the  value  of  the  reflection  coefficient  averaged  over  the  illuminated 
hemispherical  surface: 


<R> 


d(cos20)  R  = 


\  d  ( cos2  T )  R  , 

Jo 


(2.23) 


where  6  is  the  angle  between  the  surface  normal  and  the  wave  vector  of  the  inci- 

1  -  ,2  -  .2 

dent  radiation,  T  =  -s-  !T  -  0 ,  and  R  =  r  +  r  ,  with  r  and  r  the 
1  p  n  p  n 

Fresnel  reflection  coefficients  for  the  two  polarizations.  By  numerical  evaluation 
of  (2.23)  for  specific  values  of  n,  it  is  not  difficult  to  show  that  R^  ^  (R)  in  gen¬ 
eral,  although  they  do  not  differ  greatly  for  a  number  of  different  values  of  n . 

For  example,  for  n  =  1. 27  -hi  1.37,  <R>  =  0.66  (p.  292  of  Ref.  24)  and 
=  0.72.  These  two  values  differ  by  8%,  which  is  a  typical  value  of  other 


values  of  n  considered. 
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It  is  believed  that  the  result  (2. 11)  is  correct  since  it  was  derived  from 
23,24 

the  exact  theory  of  Mie.  The  difference  between  (2.11)  and  (2.22)  is  not  im¬ 
portant  in  the  present  study  since  the  small  difference  is  overshadowed  by  other 
effects,  such  as  deviations  of  particle  shapes  and  sizes  from  the  assumed  spheres 
of  common  radius.  In  passing,  notice  that  the  phenomenological  expression  (2. 22) 
results  from  using  geometrical  optics  approximations  and  assuming  that  the  re¬ 
flections  for  both  polarizations  from  the  (curved)  differential  area  on  the  sphere 

are  the  same  as  those  of  a  plane  semi-infinite  medium  at  the  same  angle  of  inci- 

26-28 

dence.  With  these  approximations,  the  Franz,  Depperman,  Foch  surface 

waves  do  not  arise  and  surface  wave  "spray"  is  neglected.  The  exact  theory  of 
23,24 

Kie  does  contain  these  effects. 

Limiting  values  of  the  absorption  cross  section  have  been  discussed  above. 

For  example,  (2.4)  and  (2. 11)  are  the  expressions  for  a^g  for  the  two  cases  of 

lea  «  1  *ud  ka,  y0ja  »  1,  respectively.  In  some  cases  the  two  limiting 

values  are  joined  smoothly  by  a  monotonically  increasing  function  of  particle 

radius  a ,  as  illustrated  schematically  in  Fig.  1.  In  other  cases,  o^g  as  a 

function  of  a  may  exhibit  a  maximum,  as  illustrated  in  Fig.  2.  In  this  figure, 

2 

the  value  of  o^g/ira  at  ka  =  ®  was  corrected  (1  -  Rj)  =  0.75,  as  compared 
with  1  -  (  R)  =0.66  in  the  reference). 

The  value  of  o^  at  the  peak  not  only  is  greater  than  the  a  ■*  ®  limiting 

2  2 
value,  but  also  may  be  greater  than  ffa  .  In  Fig.  2,  for  example,  -  2  it  a 

at  the  peak.  An  absorption  coefficient  greater  than  the  geometrical  cross  section 
2 

17  a  is  possible  since  the  electric  field  in  the  inclusion  may  be  larger  than  the 
incident  field,  as  is  well  known  in  electromagnetic  theory  (large  fields  near  sharp 
cracks  in  dielectrics,  for  example). 


99 


n- 1.27 -1.371 


Fig.  2.  Example  of  a  cross  section  with  a  maximum  value  greater 
than  ffa*.  After  Van  der  Hulst,  Ret  24,  p.  276. 
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HI.  Absorption  Coefficients  for  Various  Types  of  Inclusions 

The  results  (2. 1),  (2. 4),  (2. 11),  and  (2.21)  will  now  be  used  to  determine 
the  effects  of  various  types  of  inclusion  in  crystals.  First  consider  the  case  of 
ka  »  1  and  fya  »  1 .  From  (2. 1)  and  (2. 11) 


P  =  tra2(l-R1)NI  , 


for  ka,  /3ja  »  1  . 


(3.1) 


This  result  for  j3  is  essentially  independent  of  temperature  in  many  cases  of  inter¬ 
est.  In  the  infrared  region,  the  frequency  dependence  of  large  metallic  inclusions 
is  either  P  ~  u  '  or  /3  independent  of  w ,  as  discussed  at  the  end  of  the  present 
section.  Near  the  reststrahl  region  of  dielectric  materials,  is  strongly  fre¬ 
quency  dependent ;  otherwise  Rj^  is  relatively  constant  in  the  infrared  region. 

If  the  absorption  coefficient  /3j  of  the  inclusion  becomes  too  large,  R^ 
approaches  unity  and  P  becomes  small.  In  the  other  limit  of  small  J3j 
(^a  «  1  «  ka ),  the  rays  traversing  the  inclusions  are  not  completely  ab¬ 
sorbed,  and  P  becomes  small  again. 

Consider  this  latter  case  of  large  inclusions  (ka  »  1)  that  absorb  only  a 
fraction  of  the  radiation  in  the  inclusion  fya  «  1 ) .  From  (2. 1)  and  (2.21), 


P  -  f  (1-Rj_)f Pi  , 


pT a  «  1  «  ka  , 


(3.2) 


where  f  =  4  tra'3  Nj/3  is  the  fractional  volume  of  inclusions. 
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Next  consider  the  small-inclusion  case  of  k  a  «  1 .  From  (2. 1)  and  (2. 4) 


9  € 


0  s 


fk 


(€  +2f+  €, 


ka  «  1  . 


(3.3) 


For  nonmetals,  where  Cj  «  typically  is  satisfied,  (2. 1)  and  (2.6)  give 


0  - 


9  Re  n 

- - y 

(€a  +  2)2  +  €j2 


for  ka  «1,  €j  £  1  ,  (3.4) 


where  the  factor  in  the  bracket  typically  has  a  value  near  unity.  For  very  strongly 
absorbing  dielectric  materials  ( co  near  the  fundamental  resonance  or  the  elec¬ 
tronic  absorption  edge,  for  example),  the  result  (3. 3)  should  be  used,  otherwise 
(3. 4)  is  more  convenient. 

The  absorption  coefficients  in  (3. 2)  and  (3. 4)  are  strongly  temperature  and 
frequency  dependent  in  general.  If  0j  is  controlled  by,  say,  the  n -phonon  sum¬ 
mation  process,  then  0  <■**  Tn  *  in  the  high  temperature  limit  and  0  decays 
exponentially  with  frequency. 

For  the  case  of  -|-  ( 1  -  )  ^  9  Re (n)  [ ( +2  )2  +  €j2  ]  2  1 ,  (3. 2)  and 

(3.4)  give 

0  -  f  0j  ,  for  both  ka  »  1  and  ka  «  1  .  (3. 5) 

For  the  cases  above  [(3.2),  (3.4),  and  (3.5)]  for  which  0  ~  0j,  a  small  fraction 
of  the  radiation  entering  the  inclusion  is  absorbed,  while  for  0  ( 1  -  )  IT  a  in 

(3.1),  essentially  all  the  radiation  entering  the  inclusion  is  absorbed,  roughly 
speaking. 
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r 


Next  consider  the  case  of  small  (ka  «  1 )  metallic  inclusions.  From 


(2. 1)  and  (2. 9), 


ft  3  2ti  . 
9  n^  co  if 


(3.6) 


This  expression  (3. 6)  shows  that  /3  increases  quadratically  with  frequency  and 
is  independent  of  temperature,  roughly. 

The  absorption  by  small  metallic  inclusions  has  a  peak  that  lies  in  the 

10 

ultraviolet  or  visible  region  typically.  The  value  of  0  at  this  peak  is  much 
greater  than  the  value  in  the  infrared.  Thus  ultraviolet  and  visible  measure¬ 
ments  can  be  used  to  verify  the  source  of  absorption  by  small  metallic  particles 
in  the  infrared. 

Consider,  as  an  example,  small  potassium  spheres  in  KC1  or  KBr. 

F  centers  can  be  transformed  to  colloidal  potassium  (small  spheres)  by  heating 
29-31 

.i  _  .  .  t  ^ A  'Ti _ _ _ £ _ _ •_ _ i _ i  t _ ..1, _ _ _  nm _ 


the  crystal. 


The  transformation  is  enhanced  by  ultraviolet  radiation.  The 


wavelengths  Aq  of  the  peaks  for  the  small  potassium  spheres  in  KC1  or  KBr 
are  0. 730  and  0. 770 |an,  respectively.^  In  the  visible  and  ultraviolet, 
r2  «  co2 ,  and  (2. 7)  and  (3. 3)  give,  with  k  =  n^  co/c 

w  2  u  2  r 

«  -  *oo  “  — +i  ■  P  [  > 

W  «H  W 


r2“2gpk 

(o>2-«02)2+  (rw02/w)2 


Sec.  C 


where 


^H+< 


9  nH3<O04  , 

C  w  2  r 

p 


Notice  that  for  fairly  narrow  lines  ( T  ^  Wq  /10 ),  P  -  P^  at  the  peak  at 
u)  =  Wg  ,  and  T  is  the  full  line  width  between  the  points  p  =  -j-  fi  ^  . 

Dividing  0^  by  the  infrared  absorption  coefficient  (3. 6)  gives 


5* __  _y 

^IR  WIR^ 


(3.8) 


For  Xg  =  2irc/o)g  =  0.75j*m,  X^  =  10.6pm,  and  T=  5x  10 ^  sec  1, 
(3. 8)  gives 


=  6.4  x  103  . 

*10.6 


(3.9) 


If  it  is  assumed  that  the  limit  of  detectability  of  color  for  the  colloids  in 
KBr  or  KC1  corresponds  to  0plc  =  0.03  cm”1,  the  value  of  0jq  ^  from  (3.9) 
is  &iq  $  =  5  X  10  ^cm  *.  Thus  the  contribution  to  p  from  potassium  colloids 
in  water  while  KC1  or  KBr  crystals  should  be  less  than  ~  5  x  10  ^cm"1,  accord¬ 
ing  to  the  order -of- magnitude  estimate.  Colored  crystals  could  have  a  greater 
contribution  to  p  from  this  source.  It  should  be  mentioned  that  impurities  in 
the  form  of  F  centers  which  give  rise  to  strong  absorption  in  the  visible  region, 
may  not  give  rise  to  detectable  absorption  in  the  infrared. 


The  usual  concept  of  reflectivity  loses  its  meaning  when  ka  «  1.  However, 
in  this  case  of  a  small  sphere,  the  incident  electric  field  is  nearly  constant  across 


the  sphere,  and  the  field  Ej  inside  the  sphere  can  be  obtained  from  simple  electro¬ 
static  theory.  The  well  known  result  is 


C-T  = 


cf 2EH 


9 


where  EH  is  the  field  in  the  host  crystal  in  the  absence  of  the  inclusion.  Thus, 
for  a  sphere  with  |  c  |  »  1 ,  the  field  inside  the  sphere  is  much  smaller  than 
that  in  the  host  material. 


Consider  the  frequency  dependence  of  1  -  Rj^  in  (3. 1).  For  the  case  of 
2  2 

r  »  co  ,  (2. 8)  gives  the  familiar  result 

‘V  _  4ira 

*  '  V*-  "  eh“  ’ 

2 

where  a  =  cop  /  4  itT  is  the  conductivity.  For  €  —  iCj  ,  the  value  of  n  is 

n=vr7Iis  /  2~ %  \  !  (1  +  i),  and  (2. 18)  gives 
’  H  ' 


With  (3. 1),  this  gives  0  ~ 

2 

For  t’ae  case  of  T  a 


1/2 

co  '  ,  as  discussed  under  (3. 1). 
TBu2  «  co2,  (2.8)  gives 


(3. 10) 


n  =  \/~P  — 

2nHcoz  rBu  / 


(3.11) 


Substituting  (3. 11)  into  (2. 18)  gives 


2nH  FBu 


r2  a-  r^2  «  w2  .  (3.  i2) 


From  (3.12)  and  (3. 1),  /3  is  independent  of  frequency,  as  mentioned  in  the  second 
case  under  (3.1).  The  temperature  dependence  of  ^  in  this  case  is  determined 
by  Fgu  iii  the  usual  case  in  which  n^  and  top  are  nearly  independent  of  T. 
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IV.  Examples 

_3 

Volume  fractions  of  various  inclusions  in  a  nonabsorbing  host  for  ^  =  10  cm 

will  be  calculated  in  this  section.  For  large  nonmetallic  particles  with  1  -  ~  1 , 

2 

(3. 1)  gives  P=  ir a  Nj.  Thus,  a  single  nonmetallic  particle  of  radius  a  =  50ilm 

3 

(with  ka,  0ja  »  1)  in  a  1cm  crystal  gives 
P  s  10  4  cm  1  . 

3  -7 

The  volume  fraction  for  this  inclusion,  f  =  (4 fra  Nj/3 ),  is  f  =  5  x  10  .  The 
-4  3 

same  value  (/3  =  10  cm  )  is  obtained  with  100  particles  per  cubic  centimeter  of 

a  =  5pm  (with,  say,  k  5  ^host  =  2  kyac  =  4ff/10pm,  so  that  ka  =  2n  »  1 

•8 

cJm(n)  -  1).  The  volume  fraction  of  inclusions  in  this  case  is  f  =  5  x  10  . 

Clearly,  such  macroscopic  inclusions  must  be  essentially  eliminated  in  order 
to  obtain  small  values  of  p.  Metallic  inclusions  (with  ka  » 1 )  are  not  as  trouble¬ 
some  as  nonmetallic.  For  a  metallic  inclusion  with  R^=  0.  95,  or  1  -  R^=  1/20, 

-4  -1 

20  times  more  inclusions  are  required  to  give  P  =  10  cm  .  Even  with  such  a 
reduction  in  absorption,  the  absorption  by  large  metallic  inclusions  is  not  negli¬ 
gible;  the  volume  fractions  in  the  two  examples  above  become  10  ^  and  10  , 
which  are  still  quite  small. 

2  2 

For  the  small  inclusion  case  of  9Re(n)[(c^+ 2)  +  ]  2:  1  in  (3. 4),  or 

for  ^(  1  -  Rj^)  ^  1  in  (3.2),  with  f  =  10  ^  and  =  10cm  *,  (3. 4)  or  (3. 2)  gives 
P  =  10  "*cni  \  For  /3j  =  lO^cm  f  =  10  Ogives  /3  =  10  Consider  the  specific 
example  of  ZnS  in  ZnSe,  which  is  of  interest  because  Kodak's  Irtran  4  is  polycrys- 
talliue  ZnSe.  Until  recently,  ZnS  was  added  to  inhibit  grain  growth.  It  is  believed 
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that  ZnS  collects  along  grain  boundaries.  With  Re  n^  =  2. 47  and  ^  =  0. 4  cm 
for  ZnS  and  n^  =  2. 89  for  ZnSe,  Ren  s:  n  a  0. 855,  and  for  ka  «  1  and  «  1 , 
(3.4)  gives 

P  2=  0.38  f  cm'1  . 

-4-1  -4 

The  volume  fraction  corresponding  to  ^  =  10  cm  is  f  =  2. 6  x  10  .  The 

results  for  fya  «  1  «  ka  from  (3.2)  are  P  =  0.29  f  cm'1  and  P  =  lO'^cm"1 
-4 

for  f  =  3. 4  x  10  .  Thus,  <*<0. 01  to  ~0. 03  volume  percent  of  ZnS  in  ZnSe  gives 
P  =2  10"4  cm'1. 

The  present  theory  indicates  that  the  10. 6/im  value  of  P  ss  0. 1  cm’1  for 
Irtran  4  is  not  a  result  of  ZnS  in  the  crystal.  Since  /3.  =  0.4  cm1,  approximately 
25  percent  ZnS  would  be  required  to  give  P  =  0. 1  cm"1.  Even  the  doped  samples 
did  not  contain  such  a  large  fraction  of  ZnS. 

Next  consider  the  case  of  very  small  metallic  inclusions.  For  w  £r  2  x  1014 
rad/sec,  T  as  Tgu  —  5  x  1014sec  1  (for  a  =  20 A),  =  5x  lO1^  rad/sec, 

k  —  2  ff/  10. 6fim,  and  €„  =  4,  (3. 6)  gives 

VoC  fl 

P  as  2000  f  cm-1  . 

-8  -4  -1 

For  a  very  small  fractional  volume  of  f  =  5  X  10  ,  this  gives  P  =  10  cm 

for  the  small  metallic  inclusions. 

As  a  final  example,  consider  the  problem  of  a  surface  contaminated  by, 
say,  a  polishing  compound.  Since  polishing  compounds  contain  light  elements 
such  as  oxygen,  aluminum,  and  carbon  (in  order  to  obtain  hardness),  their  funda¬ 
mental  resonance  frequencies  are  high,  that  is,  near  the  10. 6 Jim  frequency. 
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Thus,  their  absorption  coefficients  are  large  at  10. 6 Jim.  For  example, 

extrapolating  the  A^O^  data  on  the  range  4 Jim  <  X  <  5. 5  Jim  to  10. 6jim 
2  4-1 

gives  j3j  =  10  -  10  cm  .  Even  though  the  extrapolation  procedure  gives  a 
highly  uncertain  value  of  /3j ,  a  value  of  /3j  =  100  cm  *  seems  quite  reasonable; 
the  values  of  J3  for  LiF,  MgO,  and  GaP  at  10. 6 Jim  are  50,  60,  and  123  cm 
respectively. 


As  a  simple  specific  model,  consider  particles  of  radius  a  =  ji  m  in  a 

surface  layer  of  thickness  LSU  "  2a*  The  fraction  volume  in  the  surface  layer 

-2  -1 
is  taken  as  fSu  =  10  ,  and  the  value  of  /3j  is  assumed  to  be  100  cm  .  The 

value  of  j3  for  the  surface  layer,  from  (3. 5),  is  )3gu  2:  fgi;/3j  =  1 .  The  fraction 

of  the  energy  absorbed  in  the  surface  layer  is  /3gu  Lgu  =  5  x  10 Thus,  the 

-4 

fractional  absorption  at  two  surfaces  is  ~10  . 


2 

In  spite  of  the  fact  that  the  cross  sections  are  often  much  less  than  ffa  , 
they  are  quite  large  in  a  number  of  important  cases.  As  a  point  of  reference, 
recall  that  the  real  and  imaginary  parts  of  the  wavevector  are  equal  for 

J3j  =  4tt/10.  6jum  =  10* cm  *,  roughly,  at  10. 6jim.  See  Eq.  (2.5).  For 

23  -3  -19  2 

N  =  10  cm  ,  the  corresponding  cross  section  is  a  =  j3  /  N  =  10  cm  . 

-19  2  4  -1 

These  large  values  of  a  =  10  cm  or  /3  =  10  cm  are  typical  of  values 
at  the  fundamental  resonance  in  ionic  crystals. 


-4  -1  -7  -8 

For  a  value  of  0  =  10  cm  and  a  volume  of  f  ~  10  to  10  ,  which 

corresponds  to  ~10^  to  10^  impurity  atoms  per  cm^,  the  value  of  the  cross 

section  per  impurity  atom  is  ~10  */10^  to  10  */10^  =  10-1^  to  10  2^cm2. 

Thus,  the  average  cross  section  per  atom  in  a  highly  absorbing  inclusion  is 

large,  approaching  that  of  an  atom  in  a  solid  absorbing  at  the  fundamental 


resonance. 


V.  Material  Failure  From  Local  Heating 


The  heating  of  macroscopic  inclusions  can  give  rise  to  localized  regions 
of  high  temperature  that  can  cause  material  failure  when  the  intensity  i  •  large. 
This  is  not  a  nonlinear  effect,  but  is  usually  important  only  at  relatively  high 
intensities.  Consider  the  same  model,  used  in  the  previous  sections,  of  a 
uniformly  heated  spherical  inclusion  of  radius  a  in  an  infinite  host  crystal 
that  does  not  absorb  the  laser  radiation.  For  times  short  with  respect  to  a 
characteristic  time  for  heat  to  diffuse  by  a  distance  a ,  very  roughly  speaking, 
most  of  the  energy  absorbed  by  the  inclusion  remains  in  the  inclusion.  Thus 
the  temperature  in  the  inclusion  increases  linearly  with  time.  For  times  large 
with  respect  to  the  characteristic  time,  part  of  the  absorbed  energy  has  dif¬ 
fused  into  the  host  crystal,  and  the  temperature  rise  in  the  inclusion  is  con¬ 
siderably  less  than  the  value  obtained  by  neglecting  diffusion.  Thus,  a  given 
amount  of  energy  will  cause  a  greater  temperature  rise  if  it  is  deposited  in  a 
time  that  is  short  with  respect  to  the  c'uaracteristic  time  than  if  deposited  in 
a  time  long  with  respect  to  the  characteristic  time;  the  local  heating  of  macro¬ 
scopic  inclusions  is  a  more  severe  problem  in  high-intensity  short-nulse  systems 
than  in  low-intensity  long-pulse  or  cw  systems  of  equal  average  intensity. 

The  criterion  for  failure  of  the  window  material  depends  on  the  details,  of 
the  laser  system  and  the  type  of  material  and  inclusions.  Since  there  is  no  uni¬ 
versal  criterion,  it  will  be  assumed  that  a  temperature  rise  of  1,000  K  consti¬ 
tutes  failure.  This  is  a  reasonable  choice  for  the  following  reasons:  This 
temperature  is  approximately  the  correct  value  for  melting  temperatures  and 
fracture-inducing  temperatures.  The  latter  have  typical  values  of  the  order 
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of  32  *34  (j^/aE,  where  <rf  is  the  material  strength,  a  the  linear  thermal 
expansion  coefficient,  and  E  the  Youngs  modulus.  For  =  10^  psi,  a  =  10  ** 
and  E  =  10  psi,  the  temperature  corresponding  to  fracture  is  1,000 K.  Heats 
of  fusion  have  typical  values  corresponding  to  several  hundred  degrees  Kelvin. 
Ac  1, 000  K  above  ambient  temperature  in  materials  which  do  not  melt,  the 
ionic  diffusion  may  be  important.  Although  order -of- magnitude  accuracy  in 
temperature  usually  is  not  sufficient,  order -of- magnitude  accuracy  of  intensi¬ 
ties  corresponding  to  failure  is  all  that  can  be  expected  at  present,  and  this  is 
often  sufficient.  Since  the  present  interest  is  in  this  failure  intensity  If  and 
the  failure  temperature  is  linearly  related  to  1^,  the  value  =  1,000K 
should  be  sufficient  for  present  purposes. 

Some  differences  between  the  failure  criteria  for  high-power  windows 
and  for  previous  studies  of  glasses  for  glass -host  lasers  (at  1.06/im,  for 
example)  should  be  mentioned.  In  the  latter,  highly  focused,  single  nano¬ 
second  pulses  are  used,  and  breakdown  occurs  in  a  small  volume.  In  the 
former,  all  of  the  window  is  illuminated,  apart  from  edges,  .  v  'rt- structure 
shadows,  Cassegranian-optics  "shadows",  etc. ,  and  operati  ations  range 
from  microsecond  pulses  with  repetition  rates  of  a  few  hundred  pulses  per 
second,  to  ~  1  sec  pulses  with  several  seconds  between  pulses  in  order  to 
allow  the  window  to  be  cooled,  to  continuous  operation. 

Thus,  two  important  features  of  high-power  laser-window  failure  are 
that  failure  of  the  weakest  part  of  the  window  can  constitute  system  failure, 
and  that  fatigue  and  other  multiple-pulse  effects  must  be  considered  since  re¬ 
peated  pulses  must  be  withstood.  As  an  example,  in  a  single -pulse  measure¬ 
ment,  a  laser  glass  conceivably  could  melt  locally  and  recrystallize  without 
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leaving  detectable  damage.  For  a  window  in  a  pulse -operated  system,  the 
local  absorption  coefficient  could  be  changed  by  the  higfr  temperature  associ¬ 
ated  with  the  first  pulse  or  the  first  1:  pulses,  thus  causing  increased  absorp¬ 
tion  in  subsequent  pulses  with  eventual  failure. 

Bloembergen  has  suggested  that  local  field  enhancement,  such  as  that 
occurring  at  the  edge  of  a  crack  in  a  material,  may  give  rise  to  local  intensi¬ 
ties  up  to  100  times  greater  than  the  nominal  external  intensity.  This  effect 
could  lower  the  values  of  the  breakdown  intensities  calculated  below  by  a  factor 
of  >■*'100. 

The  temperature  rise  for  T  <  T^ ,  where  Tm  is  the  melting  tempera¬ 
ture,  can  be  calculated  simply  for  the  following  model.  The  spherical  inclusion 
of  radius  a  has  temperature -independent  values  Cj  and  Kj  of  heat  capacity  per 
unit  volume  and  thermal  conductivity.  The  host  crystal  has  temperature- 
independent  values  and  .  The  boundary  between  the  two  is  assumed  to 
be  thermally  perfect ;  that  is,  there  is  no  thermal  impedance.  Heat  absorption 
by  the  host  crystal  is  assumed  negligible. 

The  thermal  time  constant 

T  =  Ca2/4K 

is  the  time  required  for  heat  to  diffuse  a  distance  a,  roughly  speaking.  For 
short  times  t  «  Tj ,  T^,  where  the  subscripts  I  and  H  denote  the  values 
of  T  in  the  inclusion  and  in  the  host  crystal,  respectively,  the  diffusion  of 
heat  out  of  the  inclusion  is  negligible.  Thus,  the  heat-flow  equation 

-  K  V2  T  +  C  I7  =  S  ,  (5.1) 

o  t 
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in  the  inclusion,  where  I  is  the  incident  intensity  and  a  is  the  cross  section  of 
the  inclusion,  as  discussed  in  Sec.  II.  It  should  be  mentioned  that  (5. 1)  will  be 
solved  separately  in  the  inclusion  and  in  the  host  material,  and  the  appropriate 
matching  of  the  temperature  and  heat  flow  at  the  boundary  will  be  made.  Thus 
K  is  a  constant  in  each  region,  and  the  usual  spatial  derivatives  of  K  do  not 
appear  in  (5. 1). 

For  T  defined  as  zero  at  the  time  the  laser  is  turned  on  (t  =  0),  the 
solution  to  (5. 2)  is 

TI  =  3CTabsIt/4,ra3ci  »  for  1  <<:  TP  tH  *  (5.4) 

There  is  a  steady-state  solution  to  (5. 1),  which  can  be  obtained  as  follows. 
The  function 

TI  =  Tc  •  67^  r2  =  Tc  -  AT  (l>2  •  '  »  TI  •  tH  •  <5-5> 


114 


Sec.  C 


with 

AT  =  0absI/8,,aKI  ’ 


is  a  solution  to  -Kj  v  Tj  =  Sj, 


_2  2  1  3  2  dr 

V  r  =  —sr  -r—  r  - — 
2  3r  <jr 


as  is  easily  seen  by  substitution,  since 

6  . 


Here  T£  is  the  temperature  at  the  center  of  the  inclusion  and  AT  is  the  tem¬ 
perature  difference  between  the  center  and  edge  of  the  inclusion. 


The  heat  flow  across  the  surface  r  =  a  is 


Qa  =  -4ira2Kj 


(5.6) 


Physically,  this  result  indicates  that  the  heat  generated  in  the  inclusion  flows 
out  of  its  surface  in  the  steady  state. 

Next  consider  the  steady-state  solution  of  (5. 1)  in  the  host  crystal.  It  is 

well  known,  from  the  electrostatic  potential  of  a  point  change,  for  example,  that 

2  2 
v  ( 1  /r )  ~  6  ( r ) ,  where  6  is  the  delta  function.  Thus,  7  ( 1  /  r )  =  0  for  r  >  0 , 

as  is  easily  verified  by  taking  the  derivatives.  Thus,  the  solution  to  -K^v2  T^  =  0 

for  r  >  a  is 


aTa/r 


From  (5. 7),  the  value  of  Q  is 

a 


Q 


a 


4*aKHTa  . 


(5.7) 


(5.  8) 
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where  Tfl  is  the  temperature  at  r  =  a .  Equating  the  values  of  in  (5. 6)  and 
(5. 8)  and  using  (5. 3)  gives 


Ta  *  °abSI',4,'aKH  • 


From  (5. 5)  with  x"  =  &  " 


T  =  T  -  AT  =  T 
a  c  c 


^abs* 
8ffa  ICT 


(5.9) 


(5. 10) 


From  (5. 9)  and  (5. 10),  the  steady  state  value  of  T£  is 


T  =  AT 
c® 


r  2KI  \ 

3<W 

1  _ 

1 

1 2 

,  1  \ 

irz  ~ 

eff 

T 

+Kr) 

The  first  term  in  (5. 11)  is  the  temperature  rise  in  the  inclusion  and  the  second 
term  is  the  temperature  rise  in  the  host  crystal.  The  value  of  is  determined 
by  the  smaller  of  and  Kj  when  one  K  is  much  larger  than  the  other.  Thus 

3 

for  a  metallic  inclusion  in  an  infrared  window  material. 

The  results  (5.4)  and  (5. 11)  are  illustrated  schematically  in  Fig.  3  for  the 

three  cases  of  Kj  «  j  ,  Kj  =  ,  and  ^  Kpj  «  Kj  in  the  (a),  (b),  and  (c) 

parts  of  the  figure,  respectively.  For  all  three  cases,  Tc  increases  linearly 

with  time  according  to  (5.4)  for  t  «  T  Consider  the  values  of  T,  for 

t  «  T  For  Kpj  -*  ®,  (5.9)  gives  T&  =  0 ;  the  infinitely  conducting  host 

maintains  the  temperature  of  the  surface  of  the  inclusion  at  T  =  0 .  In  this  case, 

2 

(5. 11)  and  (5.5)  give  Tca>  =  AT  =  -j-T^,  and  the  T^  asymptote  lies  below 
T  j  in  Fig.  3a.  In  Fig.  3b,  for  Kj  =  K^,  Tc<b  =  2  T^  . 


/ 
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Fig,  3.  Schematic  illustration  of  the  temperature  at  the  center  of  a  uniformly  heated 
spherical  inclusion  for  three  limiting  values  of  K„/K..  See  Eqs.  (5. 4)  and  (5. 11). 
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For  Kh  -*  0,  (5. 11)  gives  Tcoo  00  .  There  is  no  finite  steady-state 
solution  since  the  host  crystal  conducts  no  heat  out  of  the  inclusion.  The  tem¬ 
perature  of  the  inclusion  rises  linearly  with  time  for  all  times  according  to 
(5.4).  For  KH  large  but  finite,  T^  is  large  as  illustrated  in  Fig.  3c. 

Surface  heating  of  inclusion.  Next  consider  the  case  in  which  the  heat  is 

generated  near  the  surface  of  the  inclusion,  rather  than  uniformly  throughout 

its  volume.  Such  is  the  case  for  typical  metallic  inclusions  in  dielectric  hosts 

or  in  the  extreme  case  of  /3  a  »  1  for  other  types  of  inclusions  in  dielectric 

4  -1 

hosts.  As  an  example  of  the  latter,  for  0  =  10  cm  and  a  =  10  Jim,  /3a  =  10. 

_  35 

The  skin  depth  0  for  a  metal  is 

-1/2 

6  =  c  (27 two)  (5. 12) 

in  Gaussian  units,  where  c  is  the  velocity  of  light,  to  is  the  angular  frequency, 

o  10 

and  a  is  the  conductivity.  A  typical  value  for  6  is  40  A  for  Cu  with  A  =  10  |im. 

Energy  is  absorbed  in  the  metal  at  a  rate  which  decreases  exponentially  with  dis- 

2 

tance  into  the  metal  [S  «  E  «  exp(-2x/6)]. 

For  spherical  inclusions  with  radius  a  »  6 ,  we  assume  that  heat  is  gen¬ 
erated  uniformly  within  a  layer  of  thickness  6  over  the  entire  surface  of  the  in¬ 
clusion.  This  is  a  good  approximation  for  a  «  A  because  the  electric  field  is 
nearly  constant  over  distances  of  order  a .  For  larger  radii,  there  will  be  local 
hot  spots  over  the  surface  not  only  on  the  front  surface  but  also  on  the  rear  sur¬ 
face  , 
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There  are  three  characteristic  times  of  interest.  First, 


T6  =  4Cj  62/ wKj  (5.13) 

is  roughly  the  time  in  which  heat  diffuses  out  of  the  skin  depth,  assuming  negli¬ 
gible  diffusion  into  the  host  for  small  time  since  (dielectric)  «  Kj  (metal). 
Second, 

Ta  =  4CIa2//9irKI  (5.14) 

is  the  time  in  which  heat  diffuses  from  the  inclusion  surface  to  the  center, 
roughly  speaking.  Third, 

tH  »  CHaV3KH  (5.15) 

is  roughly  the  time  in  which  heat  diffuses  a  distance  equal  to  the  radius  of  the 
inclusion  into  the  host.  The  values  of  the  numerical  coefficients  in  (5. 13)-(5. 15) 
are  chosen  for  later  convenience. 

The  temperature  at  the  inclusion  surface  will  now  be  determined  for  four 
time  scales:  For  t  «  Tg ,  heat  does  not  diffuse  out  of  the  skin  depth,  and  the 
surface  temperature  is  linear  in  time.  At  later  times,  but  with  t  «  r  ,  heat 

is  diffusing  out  of  the  skin  depth  and  the  thickness  of  the  heated  region  increases 

1/2  1/2 
as  t  .  Thus,  the  surface  temperature  increases  as  t  .  For  still  larger 

times  Ta  «  t  «  Tj_j,  the  temperature  within  the  sphere  reaches  a  spatially 
constant  value  and  increases  linearly  with  time.  Finally,  for  t  ^  ,  heat 

diffusion  into  the  host  becomes  important  and  the  temperature  within  the  inclu¬ 
sion  reaches  an  equilibrium  value. 
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For  t  «  Tg  and  6  «  a ,  the  spherical  shape  of  the  inclusion  surface 

is  not  important,  and  we  can  solve  the  easier  problem  of  heat  generation  in  a 
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thin  plane  slab.  The  Laplace  transform  of  the  diffusion  equation  (5. 1)  is 


d  T  ( *».  P)  -  q^T 
dx2 


(5. 16) 


In  the  present  problem,  heat  is  generated  at  a  rate  S  =  I  CTa^s  /4ira  6  in  the 
slab  J  <  x  <  6 ,  where  x  =  6  corresponds  to  inclusion  surface  r  =  a .  Equa¬ 
tion  (5. 16)  inside,  outside,  and  within  the  skin  depth  of  the  inclusion  becomes 


d^  Tj  (p)  2  /D) 

- \ - qI  TIP  "  0  *  x<0  (5.17a) 

dx 


d2TH(p) 

dx2 


-q  2  T  ^  =  0 
qH  lH  U 


>  6 


(5. 17b) 


d2  Tc 


dx 


(P)  2 

1 - qI 


X  (p)  -  ^ 

6  KjP 


0  <  x  <  6 


(5. 17c) 


Equation  (5. 17)  has  the  solution 


TI<P>  =  AeV  • 

x  <  0 

(5. 18a) 

/  x  "qHX 

TH(p)=  Be 

x  >  6 

(5. 18b) 

s  qTx  -qTx 

Tc  (p)  =-f-  +  Ce*  +  De  1  , 

0  <  x  <  6  . 

(5.18c) 

P  cT 
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The  constants  A  through  D  are  determined  from  the  boundary  conditions  that 

J  'T' 

T  and  K  -~r—  are  continuous  at  x  =  0  and  x  =  a : 
dx 


A  =  -J-  +  C  +  D 

P  Cj 

(5. 19a) 

qjA  =  qjC  -  qjD 

(5. 19b) 

n  "V  S  ,  p  <1ia  ,  n  "qia 

Be  —  — * —  +  C  e  +  De. 

p2Cj 

(5. 19c) 

"Qua  Qt®  "Qra 

-KHqHBe  ”  =  K^Ce  -  K^De  1  . 

(5. 19d) 

The  temperature  at  the  inclusion  surface  is  given  by  T(p)  =  B  exp (-  qH$ ). 
Solving  (5. 19)  for  B  gives 


T(x=6,p) 


4.  /vi  X- 


Inverting  and  keeping  the  dominant  term  for  t  small  gives 


T(x=6,t) 


_S_  /ciki  ( 

CI  ) 


tCCTj  . 


Notice  that  for  CjKj  »  (appropriate  for  metal  in  dielectric) 

T(x=6,t)  =  St/Cj  as  one  would  expect.  From  this  last  equation  with 
2 

S  =  I  oaj3a/4tra  6,  the  surface  temperature  of  the  inclusion  for  small 
times  is 


Ta  =  C7 


la 


abs 

T 


I  (^Kj  +  ZC^K^)  4ira26 


t  ^ 


1  gabs 
4ira26C¥ 


t«  Te 


(5.20) 


In  (5. 20)  and  In  equations  below,  the  approximate  equalities  are  valid  for 
CjKj  »  C^Kpj.  For  Tg  «  T  «  T&,  we  can  still  use  the  approximation  of 
heat  generation  in  a  plane  slab  with  the  heat  generating  region  approximated  by 
a  delta-function  source  at  x  =  0 .  Solutions  to  the  heat -diffusion  equation  (5. 16) 
are  now 

x  <  0  ,  (5.21a) 

x  >  0  .  (5.21b) 


Tj(p,x)  =  Ae 


V 


TH(p,x)  =  Be 


V 


The  boundary  conditions  are  that  Tj  =  at  x  =  0  and  that  the  heat  flux  away 

2 

from  x  =  0  is  lo^/ 4nai  : 


i. 

p  4ira2 


(5.22) 


Solving  (5.22)  for  A  or  B  gives 


T(x=0,p) 


la 


abs 

T 


4ira 


1 

“375"  * 

P 


(5.23) 


The  inverse  transform  of  (5. 23)  is  the  temperature  at  the  inclusion  surface 


T 

a 


_ 2_ 

S7  (/cJJk^ 


I(Tabs  tl/2 
4wa2 


Tg  «  t  «  Ta  .  (5.24) 


When  t  »  t  ,  the  temperature  inside  the  spherical  inclusion  reaches 

a 

a  spatially  uniform  value.  The  present  problem  then  is  equivalent  to  uniform 
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heat  generation  within  an  inclusion  of  infinite  conductivity,  and  we  can  use  the 
previously  calculated  result  (5. 4)  to  obtain 


T  = 


31  Pahs* 
4ffa3CT 


Ta  «  t  CC  th  .  (5.25) 


For  t  »  ,  the  temperature  inside  the  spherical  inclusion  reaches  an 

equilibrium  value,  and  the  temperature  in  the  host  is  given  by  Tu  =  T  ( a/r ). 

H  Si  co 

The  constant,  ,  is  determined  by  equating  the  heat  flux  away  from  the 
inclusion  to  I  o  .  , 


Iaabs=  4lTa  KH  Ta J*  ’ 


/a 


or 


T 

a  oo 


la 


abs 


4traK 


H 


'  »  TH  . 


(5.26) 


This  result  is  equivalent  to  (5. 11)  with  K.  =  ®  (so  that  the  temperatures  at  the 
edge  and  center  of  the  inclusion  are  equal). 

The  temperature  at  the  surface  of  the  spherical  inclusion  is  sketched  in 

Fig.  4  for  the  various  time  regimes.  The  extrapolated  low-temperature  linear 

1/2 

time  dependence  of  the  temperature  intersects  the  t  '  curve  at  t  =  4 . 

Equating  (5. 20)  and  (5. 24)  gives 

4  =  T6  a  4Cj  6 2  /  it  Kj  .  (5.27) 
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Fig.  4.  Schematic  illustration  of  the  temperature  at  the  surface  of  spherical 
inclusion  with  uniform  heat  generation  within  skin  depth  6  «  a  at  the  surface. 
The  approximation  CjKj  »  is  made  in  the  values  of  the  parameters.  See 

Eqs.  (5.26W5.31). 
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At  t  =  tj ,  these  extrapolations  give 


Ti  = 


16a.  , 

abs  1 

Tr  TC 


1 6  o 


ltZaZ  RI  (/CpC[+^CHKH  )  ir^a^K, 


abs 


(5.28) 


1/2 

At  t  =  t2 ,  the  extrapolated  t  curve  intersects  the  next  linear  region  given 
by  (5. 25).  Equating  the  temperatures  in  (5. 24)  and  (5. 25)  at  ,  we  find 


*2  = 


a  Jr2 
4  a  C| 


9ff(/C^K^+v/Cp^r 


(5.29) 


At  t  =  t2  ,  (5.29)  and  (5.25)  give 


T2  ~  73 


I  o  .  C. 
abs  I 


Io 


3 ir^ a  (/ CH Kh  +  / Cj Kj )^  3ir*aKj 


abs 


(5.30) 


At  t  =  tg ,  the  extrapolated  linear  curve  (5. 25)  crosses  the  equilibrium  value 
(5. 26).  Equating  (5. 25)  and  (5. 26)  yields 


Cj  a 

l3  "  3K~ 


and 


H 


T,  =  tq 

J  3.00 


-1  T 

cH  H 


(5.31) 


»-7 


For  absorption  in  a  metallic  inclusion  at  X  =  10  Jim,  with  6  =  4x  10  cm, 
2J/cm3K,  a  =  lOfim,  K 
(5.27),  (5.29),  and  (5.31)  give 


Cj  =  2J/cm3K,  a  =  lOym,  Kj  =  2W/cmK,  and  KH  =  10'1  to  10"3W/cmK, 
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tj  2i  2  x  10  ^  sec  ,  tj  —  10  7  sec 

tg  —  10  ^  to  10  ”  sec. 


(5.32) 


The  values  of  tj  and  are  shorter  than  pulse  lengths  of  current  high -power 
pulsed  CO2  lasers.  For  a  pulse  length  of  ljisec ,  the  temperature  is  given 
by  (5.25). 

The  various  time  regimes  in  Fig.  4  depend  on  inclusion  size  in  such  a 

way  that  if  the  pulse  length  is  fixed,  long  times  in  Fig.  4  are  associated  with 

small  inclusions.  For  example,  the  time  t^  at  which  the  temperature  reaches 

1/2 

equilibrium  is  equivalent  to  an  inclusion  of  radius  a  =  (3K^t/C.)  for  a 
pulse  of  duration  t . 

The  above  results  also  apply  to  dielectric  inclusions  with  a  »  1 ,  as 

well  as  to  the  metallic  inclusion  used  in  the  examples.  However,  in  the  case  of 

strongly  absorbing  dielectrics  the  following  differences  should  be  noted:  First, 

the  approximation  Kj  »  kh  is  not  generally  valid  and  Kj  2:  is  perhaps 

more  generally  true.  Second,  the  radiation  penetrates  to  a  greater  depth  into 

4  -1 

the  inclusion.  Generally  /3  <  10  cm  so  that  the  radiation  penetration  depth 
is  usually  greater  than  lfim,  as  compared  to  ~40  k  for  good  conductors.  As  a 
result,  the  first  linear  region  of  Fig.  4 ,  (5. 20),  which  was  only  of  academic 
interest  in  the  case  of  metals,  will  be  important ,  and  the  size  at  which  volume 
heating  rather  than  surface  heating  applies  will  occur  at  a  larger  radius. 

For  repetitively  pulsed  lasers,  the  above  results  may  be  applied  directly 
to  each  pulse  if  the  temperature  relaxes  to  the  initial  value  between  pulses.  In 
the  case  of  a  fast  repetition  rate,  the  temperature  will  not  relax  and  the  failure 
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temperatures  will  be  reached  at  a  lower  intensity  than  calculated  for  single 
pulse.  This  will  occur  for  two  cases :  First,  if  the  time  between  pulses  Is 
less  than  th  ,  heat  will  not  diffuse  away  from  the  inclusion  before  the  next 
pulse.  Second,  for  a  large  density  of  inclusions,  heat  that  diffuses  away 
from  one  inclusion  appears  at  another,  and  the  background  temperature  of 
the  entire  window  is  increased.  This  occurs  when  the  time  between  pulses 
is  greater  than  about  / 4K^,  where  l  is  the  distance  between  inclu¬ 
sions,  but  less  than  the  time  constant  for  cooling  the  window. 


In  the  following  contractual  period  the  dielectric  or  metallic  inclusion 
temperature  will  be  calculated  as  a  function  of  inclusion  radius  for  various 


pulse  durations.  The  heating  of  platinum  inclusions  in  a  glass  host  by  30  nsec 
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duration  pulses  has  also  been  considered  by  Hopper  and  Uhlmann.  They 


assumed  rj^  =  ira2  €x(«x 


emissivity)  which  is  not  valid  for  small  in¬ 


clusions  (long  times).  Using  the  absorption  cross  sections  discussed  in 


Sec.  II,  our  preliminary  results  indicate  that  the  plot  of  temperature  versus 
size  in  their  example  does  not  have  a  maximum  as  they  found,  but  instead  T 


is  a  monotonically  increasing  function  of  a . 


Preliminary  calculations  show  that  micron-size  metal  or  dielectric  inclu¬ 
sions  can  typically  be  heated  to  the  failure  temperature  of  1,000  K  by  a  single 

2  7  2 

microsecond  pulse  with  It  less  than  50J/cm  (1=5x10  W/cm  with  t  =  lj*sec). 

In  cases  of  optimum  heating,  local  material  failure  can  occur  at  microsecond 

2 

duration  -  pulse  energies  is  low  as  5J/;;r-:  .  Even  though  the  failure  tempera¬ 
ture  is  reached  in  repetitively  pulsed  systems  of  these  energies,  cw  systems 
with  the  same  average  intensity  will  only  heat  the  inclusions  to  about  1 K . 


;-r- 
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As  mentioned  by  Bloembergen,  local  electric  field  enhancement 
for  inclusions  near  imperfections  such  as  cracks  in  the  host  material  could 
lower  the  calculated  material  damage  pulse  energy  by  as  much  as  a  factor 
of  100. 
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VI.  Conclusions 


Small  amounts  of  macroscopic  inclusions  in  or  on  the  surface  of 

transparent  crystals  can  cause  substantial  absorption.  For  example,  a 

“7  “8  15 

fractional  volt  me  of  f  ~  10  to  10  ,  which  corresponds  roughly  to  ~10 

to  1016  impurity  atoms  per  cm'*,  can  give  rise  to  0  =  10  *cm  1  for  some 

types  of  inclusions.  The  corresponding  effective  average  cross  sections 
-19  -20  2 

per  atom,  ~10  to  10  cm  ,  are  quite  large. 

Simple  expressions  for  the  absorption  cross  sections  of  spherical  in¬ 
clusions  were  derived  for  the  cases  of  ka  «  1  [  Eq.  (2.4)]  and  of  ka  »  1 
[  Eq.  (2. 16)] .  The  results  were  used  to  derive  expressions  for  0  for  the 
four  cases  of  large  inclusions  ( ka  »  1 )  of  weak  and  strong  absorbers  and 
of  small  inclusions  ( ka  «  1 )  of  dielectric  and  metallic  particles.  The 
results  are  summarized  in  Table  I,  which  contains  expressions  for  0,  typi¬ 
cal  temperature  and  frequency  dependences,  and  typical  numerical  results. 

The  result  that  the  cross  section  (2.6)  for  ka,  0^a  »  1  contains  , 
rather  than  (R  ) ,  is  of  general  interest  in  the  theory  of  scattering. 

The  local  heating  of  macroscopic  inclusions  is  a  far  greater  problem  in 
high-intensity  short-pulse  systems  than  in  low-intensity  long-pulse  or  cw  sys¬ 
tems  having  equal  average  intensity.  Preliminary  results  indicate  that  material 

2 

failure  can  occur  for  pulse  energies  as  low  as  5J/cm  with  microsecond  pulse 
durations.  Bloembergen  has  suggested  that  local  electric  field  enhancement 
near  imperfections  such  as  cracks  in  the  host  material  could  lower  the  calcu¬ 
lated  material  damage  pulse  energy  by  as  much  as  a  factor  of  100.  The  deter¬ 
mination  of  the  thresholds  for  damage  in  typical  window  materials  will  be  further 
explored  during  the  next  contractual  period. 


Table  I.  Tabulation  of  expressions  for  /3 ,  typical  T  and  a:  dependence  of  /3 ,  and  typical  values  of 
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D.  COLLECTION  OF  EXPERIMENTAL  RESULTS  FOR  0  (to) 


Information  on  the  frequency  dependence  of  the  absorption  coefficient  has 
been  collected  from  the  literature.  All  these  results  have  been  plotted  on  the 
same  scale  for  convenience  of  comparison.  Figs.  1-17  contain  all  the  results 
collected  to  date.  Figs.  18-22  give  comparisons  of  alkali  halides,  alkaline 
earth  fluorides,  zincblend  -  structure  semiconductors,  and  type  IV  elemental 
semiconductors  plus  BaTiO^,  respectively.  Fig.  23  shows  a  comparison  of 
one  material  from  each  of  the  above  groups,  plus  Ba  TiO^ .  In  this  last  curve, 
the  elements  in  KBr,  ZnSe,  and  Ge  are  from  the  same  row  in  the  periodic 
table.  The  lozenges  denote  the  10. 6  Jim  values  from  references  39  and  43. 


For  Figs.  1-7,  the  solid-curve  data  was  presented  in  the  original  re'^r 
ences  as  graphs  of  conductivity  2a  vs  a)  or  imaginary  part  of  the  dielectric 
constant  Cj  vs  to.  According  to  the  relations. 


o  =  217(20) 

•  n 


t 


the  values  of  the  index  of  refraction  nr  (to)  are  needed  to  obtain  the  values  of  0. 

Since  the  values  of  n  (to)  used  in  these  references  were  not  given,  the  following 

values  (lowest-wavelength  values  from  AIP  Handbook)  were  used  formally: 

n  =1.5  for  KBr;  n  =1.2  for  KC1  and  NaCl;  n  =1.1  for  UF ;  n  =1.4  for 
r  r  r  r 

BaF2  J  nr  =  1*4  for  SrF2  ;  and  nr  =  1. 3  for  CaF2  .  Notice  that  this  solid  -curve 
data  is  in  the  vacinity  of  the  main  resonance,  where  nr  is  a  strong  function  of  to. 
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2  3  4  5  6  7  8 

W/COf 

Fig,  2.  Experimental  frequency  dependence  of  the  infrared  absorption 
coefficient  0  for  KC1  after  Horrigan  and  Deutsch  (+,  x,®  .O )  Ref,  39, 

Smart  et  al  (-  •- )  Ref.  41 ,  and  Genxel  (— )  Ref,  40. 
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race  of  the  infrared  absorption 
Oeutsch  (+,<D )  Ret  39,  Smart  etal 
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Pig.  6.  Experimental  frequency  dependence  of  the  infrared  absorption 
coefficient  fi  for  StFj  after  Horrlgan  and  Oeutsch  (+, ©)  Ref.  39,  and 
Denham  et  al  (- )  Ref.  42, 
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Fig,  11.  Experimental  frequency  dependence  of  the  infrared  absorption 
coefficient  0  for  AlSb  after  Turner  and  Reese,  Ref.  45. 
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